BUCKLING OF LAMINATED COMPOSITE PLATES 


sy 

sudhakar sharma 




PEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOI.O0Y KANPUR 

JULY,, 1979 



BUCKLING OF LAMINATED COMPOSITE PLATES 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


By 

SUDHAKAR SHARMA 


' to the 

DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JOtY, 1979 



Pt £ " f ^ - jf^ - £ A - 6 U C: 





ii 


CERTIgICATE 

This is to certify that the v/orh BUCKLING OF LAIVIINATED 
COMPOSITE PLATES has been carried out under our supervision 
and has not been submitted elsewhere for a degree 


(N.G.RrrfflNGAR) 

Professor 

Department of Aero. Engg. 
I . I . T . Kanpur 



V 

(P.N. MURTHY ) 



Professor 

Department of Aero. Engg 
I . I .T . Kanpur 


* 



iii 


ACKlfOWIiEDGEMENyS 

It is with great' pleasure that I acknowledge my 
deep iftdebtedhass to Ecof, P.N. Murthy and Prof* NvG.R, Iyengar 
for their invaluable guidance, critical supervision and 
constant encouragement throughout the course of this work; 

I sincerely thank Dr, P,S. Kushwaha and Mr. 

M.K, Patra for several fruitful discussions. Thanks are also 
due to many of my friends for their help in the development 
and running of computer programmes in the initial stages. 

To my wife, Asha, I am indebted for her patience 
and^ willingness to make sacrifices so that this work could be 
completed. 


I anr grateful to Mr, M,K, Patra fen? his help 
during the manuBoript stage and to Mr. R.N, Srivastava for 
his fine typing- 



TABIJS OP CONma?S 





Page 

LIST OP 

TABLES 

vii 

LIST CP 

PIGHRES 

ix 

SYMBOLS 

ABB 

NOTATION 

XV 

SYNOPSIS 


xvii 

CHAPTER 

1 

INTRODUCTION 

1 - 11 

1.1 


Characterisaljion of Laminated Plates 
and literature Survey 

1 

1.2 


Buckling of Laminated Plates 

3 

1.3 


Buckling of Laminated Plates Under 
Lateral Loads 

9 

1.4 


Scope of the Thesis 

10 

CHAPTER 

2 

CHARACTERISATION OP THE LAMINATED 
COMPOSITES 

12-41 

2.1 


Introduction 

12 

2.2 


Review of Literature 

1# 

2.3 


Analysis 


2.3.1 


Pour “Layered Angle -Ply Laminate 

20 

2.3.2 

> 

Six-Layered Angle-Ply Laminate 

22 

2.4J 


Results and Discussion 

25 

2.5 


Conclusions 

27 

CHAPTER 

3 

PORMULATION OP THE PROBLEM AND METHOD 

OP SOLUTION 

42 - 65 

3.1 


Introduotion 


3.2 


Derivation) of Equilihriunr Equations 
and Natural Boundary Conditions' ^ 

' . 43 





V 


• 

Page 

3.2,1 

Assumptions 

45 

3.2,2 

S tr ain-D isplaceme nt -Curvatiir e 
Relationships 

44 

3.2.3 

Derivation of Equilitrium Equations and 
Boundary Conditions By Energy Approach 

44 

3.3 

Reduction of Three Governing Equations 
into a Single Equation 

49 

3.3.1 

Ant iayinme trio Angle-Ply Plates 

54 

3,3.2 

Antisymmetric Cross -Ply Plates 

58 

3.4’ 

Method of Solution 

61 

CHAPTER 4 

STATIC ANALYSIS OE ANTISYMMETRIC 
ANGLE-PLY PLATES 

66 - 85 

4.1 

Introduction 

66 

4.2 

Analysis 

67 

4.2,1 

Governing Equations 

67 

4.2.2 

Boundary Conditions 

68 

4.2.3 

Loading 

69 

4,2.4 

Method of Solution 

69 

4.3 

Numerical Computations 

72 

4.4 

Discussion and Conclusions 

73 

CHAPTER 5 

■ BUCKLING OE ANTISYMMETRIC CROSS- iND 
ANGLE-PLY PLATES 

86 - 146 

5.1 

Introduction 

86 

:,-:’5V2,: 

Buckling Analysis 

88 

5.2.1 

Derivation of Governing Equations 

88 

5.2.2 

Method of Solution 

89 

5.2.3 

Boundary Conditions 

90 





vi 



Page 

5.5 

Numerical Computations 

101 

5.4 

Results and Discussion 

101 

5,4.1 

Plates With Two Opposite Ed’ges Simply 
Supported and the Remaining Edges Pree 

101 

5.4.2 

Plates With Two Opposite Edges Simply 
Supported and the Remaining Edges 

Clamped 

104 

5,4.5 

Plates With Two Loaded Edges Simply 
Supported, One Fixed and the 

Remaining Free 

108 

CHAPTER 6 

BUCKLING OP ANTISIMMETRIGALLY LMIHATED, 
RECTANGULAR PLATES ON A WINKLER - 
PASOMINAK POT3NDATION 

1 47 - 168 

6«1 

Introduction 

147 

6.2 

Buckling Analysis 

148 

6.2.1 

Antisymmetrio Angle-Ply Plates 

148 

6.2.2 

Antisymmetric Cross-Ply Plates 

149 

6.3 

Results and Discussion 

150 

CHAPTER T 

BUCKLING OP ANTISIMMEGRIC CROSS- AND 
ANGLE-PLX COMPOSITE PLATES SUBJECTED 

TO VARYING UNI-AXIAL LOADS 

169 - 194 

7.1 

Introduction 

169 

7,2 

Analysis 

170 

7.2.1 

Antisymmetric Angle-Ply Laminates 

170 

7,3 

Numerical Computations 

175 

7.4 ' 

Results and Discussion 

176^ 

CHAPTER 8 

RESULTS AND CONCLUSIONS 

195 r- 201 

8.1 

General Gonclusions 

195 

8.2 

Scope for Eurther Work 

200 

REPEREHCES 


202 - 208 



LIST OF TABLES 


Yii 


Table 

Title 

Page 

2.1 

Values of and o^g foa? Various Stacking 

Sequences for Catnmon Composite Materials 

29 

4.1 

Variation of Non-dimensional Displacement 
(wEgt^/q^b^) along the Centre Line of an 
Antisymmetrically Laminated Angle -Ply Plate 
with a/b = 8.0 

76 

4,2 

Variation of Non-dimensional Moment Resultant 
(My.E 2 t/A^iqob) along the Centre-Line of an 
Antisymmetrically Laminated Plate with 
a/b =8,0 

7T 

4.5 

Variation of Non-dimensional ^Visting Moment 
Resultant (M^^E^t/A^^ q^b) along the Free Edge 
of an Antisymmetrically Laminated Plate with 
a/b =8.0 

78 

4.4 

Variation of Non-dimensional Transverse Shear 
(QyEgt /A|^q^b) along the Free Edge of an 
Antisymmetrically Laminated Angle-Ply Plate 
with a/b =8.0 

79 

5.1 

Percentage Reduction in Buckling Load Between 
Different Layers of an Antisymmetric OPRP 
Angle-Ply Square Plate as a Function of Fibre 
Orientations (Free Edges) 

115 

5.2 

Percentage Reduction in Buckling Load Between 
Different Layers of an Antisyiimetric OFRP ^^^ ^ 
Angle -Ply Square Plate as a Function of Fibre 
Orientations ( Clamped Edges) 

114 



viii 


Table 

5.3 


Title 

Percentage RecPuction in Buckling Load Between 
Different Layers of an Antisymmetric OPRP 
Angle-Ply Square Plate as a Punction of Fibre 
Orientations (Clamped -Free Edges) 


Page 


115 



ix 

LIST 05? FIGURES 

Figure Title Fage 

2,‘I Geometry of a n-Layered Laminate 31 

2.2 Geometry of a Four Layered Laminate 32 

2.3 Geometry of a Six Layered Laminate 32 

2.4 Variation of "Coupling Stiffnesses' with 
Thiclcness Ratio for a FourrDayered 

Antisymmetric Laminate 33 

2.5 Variation of “Coupling Stiffnesses with 
Thickness Ratio for a Six Layered 

Antisymmetric Laminate 54 

2.6 Variation of Coupling Stiffnesses with 
Thickness Ratio for a Six-Layered 

Antisymmetric Laminate 35 

2.7 Locus of Orthotropic Thickness Ratios 36 

2.8 Locus of Thickness Ratios for the Partial 

Elimination of Coupling Stiffness 37 

2.9 Locus of Thiclcness Ratios for the Partial 

Elimination of Coupling Stiffness 38 

2.10 Locus of Thickness Ratios for the Partial 

Elimination of Coupling Stiffness 39 

2.11 Locus of Thickness Ratios for the Partial 

Elimination of Coupling Stiffness 40 

2.12 Locus of Thickness Ratios for the Partial 

Elimination of Coupling Stiffness 41 

3.1 Co-ordinate System of the Plate 65 

4.1 Co-ordinate System.^ of the Plate 80 

4.2 Loading Arrangement 81 

4.3 Variation of Centre-Line Leflection 82 

4.4 Variation of Bending Moment Along the Centre- 

Line 83 



X 


S’igure Title Page 

4 . 5 Variation of Twisting.: Moment Along a Free 

Edge of the Plate 84 

4#6 Variation of Transverse Shear Along a Free 

Edge of the Plate 85 

5.1 Coordinates of the Plate 116 

5.2 Variation of Buckling load with Thickness 
Ratio for an ihitisymmetric Angle-Ply Square 

Plate 117 

5.3 Variation of Buckling load with Thickness 

Ratio for an Antisymmetric Square Plate 118 

5.4 Variation of Buckling load wi*th Thickness 
Ratio for an Antisymmetric Angle-Ply Square 

Plate 1 1 9 

5.5 Variation of Buckling load with Aspect Ratio 120 

5.6 Variation of Buckling load with Aspect Ratio 121 

5.7 Variation of Buckling load with Aspect Ratio 122 

5.8 Variation of Buckling load with Fibre 

Orientation 123 

5.9 Variation of Buckling load with Plate Aspect 

Ratio of an Antisymmetric Cross -Ply Plate 124 

5.10 Variation of Buckling load with Cross-Ply 

Ratio 125 

5.11 Variation of Buckling load with Fibre 

Orientation 126 

5.12 Variation of Buckling load with Fibre 

Orientation 127 

5.13 Variation of Buckling load with Fibre 

Orientation 128 

5.14 Variation of “Buckling load as a Function of 

Thickness Ratio -^ ^^^^129 



xi 

3?igure Title Page 

5.15 Variation of Buckling load as a Punction of 

Thickness Ratio 1 30 

5.16 Variation of Buckling Load v/xih Aspect Ratio 131 

5.17 Variation of Buckling Load with Aspect Ratio 132 

5.18 Variation of Buckling Load with Aspect Ratio 133 

5.19 Variation of buckling Load with Aspect Ratio 134 

5.20 Buckling load as a Function of Angle-Ply 

Orientation for a Bquare Plate Under Uniform 
Uniaxial Oompression 135 

5.21 Buckling Load as a Function of Angle-Ply 

Orientation for a Square Hate Under Uniform 
Uniaxial Compression I 36 

5.22 Buckling Load as a Function of Angle-Ply 

Orientation for a Square Plate Under Uniform 
Uniaxial Compression 137 

5.23 Buckling load as a Function of Aspect Ratio I 38 

5.24 Buckling Load as a Function of Aspect Ratio 139 

5.25 Buckling Load as a Function of Aspect Ratio I 40 

5.26 Buckling Load as a Function of Thickness Ratio 

for a Four-Layered Laminate I 4 I 

5.27 Buckling Load as a Function of Thickness Ratio 

for a Six-Layered Laminate 142 

5.28 Variation of Buckling Load as a Function of 

Cross-Ply Ratio 143 

5.29 Variation of Buckling load as a Function of 

Cross-Ply Ratio, M 144 

5.30 Variation of Buckliiig Load with Fibre 

Orientation with Number of layers and Boundary 
Conditions as Parameters 145 

5.31 Variation of Buckling load with Thickness Ratio 
with Fibre Orientation & Boundary Conditions 

as Parameters 146 



xii 

Figure Title Page 

6.11 Variation of iTondimensional Buckling Boad with 

Winkler Foundation Stiffness for a Square 
Antisymmetrio Angle-Ply Plate with Humber of 
layers (H) & Pasternak Shear Stiffness (G) as 
Parameters 166 

6.12 Variation of Nondimensional Buckling Load with 
Angle of Fibre Orientation for a Square 
Antisymmetric Plate with Number of Layers (N), 

V/inkler Stiffness (K) and Pasternak Shear 
Stiffness (G) as Parameters 167 

6.13 Variation of Nondimensional Buckling Load v/ith 

Angle of Fibre Orientation for a Square 
Antisymmetric Angle-Ply Plate with Number of 
Layers (N), Winkler Stiffness (K) and Pasternak 
Shear Stiffness (G) as Pax'ameters 168 

7.1 Geometry and Loading of the Plate 181: 

7.2 Variation of Buckling Load with Fibre 

Orientation 182 

7.3 Variation of Buckling Load with Fibre 

Orientation 163 

7.4 Variation of Buckling Load with Aspect Ratio 184 

7.5 Variation of Buckling Load with Aspect Ratio 185 

7.6 Variation of Buckling Load with Aspect Ratio 186 

7.7 Variation of Buckling Load with Aspect Ratio 187 

7.8 Variation of Buckling Load with Aspect Ratio 188 

7.9 Variation bf BuclhLing Load v/ith Thicloiess 

Ratio for a Four Layered Square Plate 189 

7.10 Variation of Buckling Load with Thickness 

Ratio for a Six Layered Plate 190 

7.11 Variation of Buckling Load with Thickness 

Ratio for a Four Layered Plate 191 



xiii 


Figure 

6.1 

6.2 

6.3 

6.4 

6.5 

6.6 

6.7 

6.8 

6.9 

6.10 


Title Page 

Geometiy of the Plate Resting on V/ihkler- 
Pastemak Foundation 1 56 

Variation ' of Buckling Load v/ith Winkler- 

Stiffness, K for a Square intisyinmetric 

Mgle-Ply Plate 157 

Variation ’ of Buckling Load v/ith Winkler 
Stiffness, K for a Square Antisymmetric 
Angle-Ply Plate 1 58 

Variation of Buckling Load with Winkler 
Stiffness, K for a Square Antisymmetric 
Angle-Ply Plate 159 

Variation of Nondimensional Buckling Load with 
V/inlcler Stiffhess for an Antisymmetric Cross- 
Ply Plate with Number of Layers (N) and 
Pasternak Shear Stiffness (G) as Paraaieters 160 


Variation of Nondimensional Buckling Load with 
Winkler Foundation Stiffness for an Antisym- 
metric Angle-Ply Square Plate v/ith Number of 
Layers (N; and Pasternak Shear Stiffness (G) 
as Parameters 

Variation of Nondimensional Buckling Load 
with Winkler Foundation Stiffness for an 
Antisymmetric Angle -Ply Square Plate with 
Fiimber of Layers (N) & Pasternak Stiffness 
(G) as Parameters 

Variation Of Nondimensional Buckling Load with 
Winkler Foundation Stiffness for' an Antisym- 
metric Angle-Ply Square Plate with Number of 
Layers, (N) and Pasternak Stiffness, (G) as 
Parameters 163 

Variation of Nondimensional Buckling Load 
with Winkler Foimdation Stiffness for a Square 
Ant isyDometric Angle-Ply Plate with Number of 
Layers (N) and Pasternak Shear Stiffness (G) 
as Parameters 164 

Variation of Nondimensional Buckling Load with 
Winkler Foundation Stiffness for a Square 
Antisymmetric Angle-Ply Plate with N-umber of 
Layers (N) & Pasternak Shear Stiffness (G) as 
■ Parameters 165 


161 


162 



S’igure 


Title 

Page 

7.12 

Variation 
Ratio for 

of Buckling load with Thickness 
a Six layered Rlate 

• 192 

7.13 

Variation 

of Buckling v/ith Cross-Ply Ratio 

193 

7.14 

Variation 

Ratio 

of Buckling load with Cross"Ply 

194 



XV 


SymOLS MB DOTATIONS 


a, b 


Dimensions of the plate 


Ai|, D.^ : Elastic constants of the laminated 

composite plate 

E^, E25 G-^2 ip • Elastic constants of an orthotropic 

lamina 


g 

k 

G 

K 


ra 


Pasternak foundation stiffness 

Winkler foundation stiffness 

gb /Dii is the non-dimensional 
Pasternak shear stiffness 

kb^/Di 1 is the non-dimensional 
Winkler foundation stiffness 

Number of buokles in the x -direction 


M 


N 


^X^ V 

N, 


X 


or 


Qi^(e) 


Ci 

R 


is the cross-ply ratio 

\ 

even 

Number of layers 
Bending and twisting moments 
In -plane normal and shear forces 
Buckling load per unit width 

Reduced stiffnesses evaluated at 
fibre orientation angle, 9 

Lateral loading 

a/b is the plate aspect ratio 

ti/t2 is the thiclaiess ratio of a 
four-layered antisymmetric laminat e 







xvi 




e 

A 

ti/'b3 is the thickness ratio of a 
six*-layered antisymmetric laminate 

Rg 



o 

o 

t2/t3 is the thickness ratio of a 
six-layered antisymmetric laminate 

•t 



« 

# 

Overall thickness of the laminate 

U, V 



p 

« 

In-plane displacements 

0 0 

U , V 


p 

0 

Mid-plane values of in-plane 
displacements 

v; 



p 

• 

Lateral displacement 

(e$, 

F° 

y’ 



Mid -plane strain components in the 
xy plane 


"y’ 

®xy^ 

« 

Strain components in the xy plane 



■^xy) 

0 

« 

Curvatures in the xy plane defined 
by Eq. (3.1b) 

e 



o 

o 

Eibce orientation angle 

0 



« 

c 

Displacement function 




0 

0 

Roots of Eq. (3.52) 

Y 



« 

0 

Roots of Eq. (5.55) 

a 



0 

« 

Numerical factca: in Eq, (7.1) 

OLvt 



9 

. • 

m7i/a 

p 



9 

9 

imi /R 



xvii 


SYNOPSIS 

SUDHAKAR SHARIO. 

Ph. D, Thesis 

Department of Aeronautical Engineering 

Indian Institute of Technology, Kanpur 

BUCKLING OP LAMINATED COMPOSITE PLATES 

The composite materials by virtue of their high 
strength-to-weight and stiffness-to-weight ratios are ideal as 
structural elements in aerospace applications. During the 
last decade much effort has been directed towards developing 
and characterising these materials. Until very recently these 
efforts have been concentrated on the characterization of the 
basic materials and the development of information on their 
behaviour. The basic difference between isotropic materials 
and fibre reinforced composites is the variety of couplings, 
for example, between in-plane loading and out-of -plane 
deflection or between applied bending load and twisting response 
that the latter exhibit. It is well-known that this coupling 
phenomenon due to lamination asymmetry can result in substantial 
reduction of effective stiffness of the plate. Conseq[uently 
some investigators have tried to eliminate the degrading effect 
of coupling by suggesting various techniques. 

Warren and Norris have shown that it is possible to 
orient the layers of multiple-layer laminates in such a way 



that the resulting elastic stretching behaviour is isotropic. 
This, however, does not remove, in general, bending and 
stretching-bending coupling. Abrams and Scheying presented 
numerous configurations for which the elastic behaviour of the 
laminate is identical to that of a homogeneous crthotropic or 
isotropic material. Bartholomew has shown that selection of 
a suitable ply stacking order enables bending, as well as in- 
plane orthotropy. A review of relevant literature reveals that 
almost no study has been made if thickness distribution of the 
laminae constituting a laminate of a given thickness can 
achieve special orthotropy. This thesis is partly an attempt 
in this direction. 

Use of laminated plates as an irapoctant structural 
element in both the aerospace and electronics industries has 
intensified the demand for the development of analytical 
methods for anisotropic structures. Of particular Interest 
are anisotropic or layered composite materials in the form 
of thin plates subjected to in-plane compressive loads. The 
primary mode of failure of these plates is buckling. 

Considerable amount of information is available on 
the buckling of orthotropic panels and on infinitely long 
anisotropic plates under special load conditions. However, 
the buckling of finite, homogeneous and layered anisotropic 
composite plates has been receiving serious attention since ‘ 
1967* Whitney and leissa presented a closed f orm solution of 



the buckling of a simply supported arbitrarily laminated 
rectangular plate subjected to biaxial compression. Their 
results indicate a strong dependence of the buckling load on 
bending-stretching coupling, Ashton and Waddoups presented 
Rayleigh-Ritz solution for buckling of arbitrarily laminated 
rectangular plates subjected to biaxial compression and in- 
plane shear loads. 

The response of arbitrarily laminated plates is 
governed by three coupled displacement equations of equilibrium. 
Most of the analyses dealing with such plates assume suitable 
functions for the three coupled displacement components. 

However, it is desirable to combine the three displacement 
equations of equilibrium into a single higher order partial 
differential equation in terms of a single function. This 
thesis is concerned with the formulation of this problem and 
its application to the buckling of laminated cross- and angle - 
ply plates under various boundary conditions. 

The problem of the buckling of laminated composite 
plates subjected to non-uniform axial loads is another area 
which has not been investigated so far. Oharais used the 
Galerkin method to determine the buckling response of simply 
supported rectangular plates under uniform normal and shear 
loads. The present study includes the buckling response of 
laminated plates under non-uniform axial load employing the 
Galerkin technique together with formulation developed in the 
present thesis. 
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The problem of the analysis of composite plates 
resting on elastic foundation is yet another area wherein 


practically very little work has been reported. Such solutions 
may find important applications in analysing the face behaviour 
of certain types of foam-filled sandwich panels and the 
structural response of wing panels resting on a cushion of air 
during flight. Turvey has presented simple exact solutions 
for uniformly loaded, simply supported antisymmetrically 
laminated cross- and angle-ply plate strips resting on V/inkler- 
Pasternak foundation. In a recent paper he has extended the 
static analysis of composite plates resting on Winkler- 
Pasternak foundation to finite cross- and angle -ply plates. 

The buckling analysis of composite plates resting on general 
elastic foundation has not been reported hitherto and the 
present thesis is partly an attempt to fill this gap. 

The thesis is divided into eight Chapters. Chapter 1 
is concerned with the review of pertinent literature on the 
characterization of laminated composites and the stability 
of rectangular laminated composite plates. 

Chapter 2 deaJ-s with characterization of laminated 
composites. Relations between thickness ratios of the laminae 
constituting a laminate of a given thickness and reduced 
stiffnesses of an crthotropic lamina are developed and conditions 
that lead to partial or complete elimination of coupling stiff- 
nesses (special crthotropy) are studied. It is shown that for 
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antisyrranetric angle-ply laminates, unique numerical values of 
the thickness ratios exist that correspond to special octhotropy 
of the laminates. These values are further shown to he indep- 
endent of material properties and fibre orientations.. The 
possibility of partial or complete elimination, of coupling 
effects for more general ts^es of laminates is also discussed. 

Chapter 5 is concerned with formulation and method of 
solution of the plate stability problems, A system of three 
coupled equations derived from energy approach and governing 
the behaviour of general unsymmetric laminates is reduced to a 
single eighth order partial differential equation in terms of 
a single displacement function. The displacement, force and 
moment resultants are also expressed in terms of the same 
function. A closed form solution is then found fcsr anti- 
symmetric cross- and angle-ply laminates. 

In the fourth chapter, a specimen problem entitled 
"Static analysis of an antisymmetrically laminated angle-ply 
plate subjected to pure bending with two opposite edges simply 
supported and the remaining two free" is solved and the results 
thereof in a special case are compared v/ith those available in 
literature to confirm the validity of the method proposed in 
Chapter 3. It is shown that the proposed method leads to a 
faster convergence. 

In the fifth chapter, buckling response of anti- 
symmetrically laminated cross- and angle-ply plates with two 
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loaded edges simply supported and various other combinations 
of boundary conditions on the remaining two edges is consi- 
dered. The effect of number of parameters namely thickness 
ratio of laminae constituting a laminate of a given thickness, 
cross-ply ratios, number of layers and fibre orientations etc. 
on the buckling response of plates is discussed. 

Chapter 6 is concerned with the buckling response 
of antisjrrametrically laminated cross- and angle -ply plates 
with two opposite loaded edges simply supported and various 
other boundary conditions on the remaining edges resting on 
general elastic (Winlcler-Pasternak) foundation. Typical 
results for a number of different plate lay-ups and Winkler- 
Pasternak values are presented in dimensionless graphical form 
The effect of material coupling and foundation stiffnesses on 
the buckling response of plates is discussed. It is shown 
that material coupling and foundation stiffnesses can affect 
the buckling response of plates significantly. 

In the seventh chapter, buckling response of simply 
supported antisymmetric cross- and angle -ply plates subjected 
to pure in-plane bending is discussed. The problem is solved 
using Galerkin method together with the formulation developed 
in Chapter 3. The results in a special case (with uniform uni 
axial loads) are c crap ared with those available in literature. 

The thesis concludes with summing up of the majcsr 
conclusions drawn from these studies together with suggestions 
for further work. ' 



CHAPTER 1 


INTRODUCTION 

1.1 CHARACTERISATION OF LAMINATED PLATES MD LITERATURE SURVEY 

Daminated plates constitute an important structural 
element in both the aerospace and electronics industries. In 
aerospace applications where weight savings are of paramount 
importance, the advent of fibre-reinforced composite materials 
such as carbon/epoxy and graphite/epoxy has resulted in a 
dynamic increase in the use of laminated fibre -reinforced plates 
and other structural shapes. The composite materials are 
typically a combination of a usually light, wealc and flexible 
binding material, v/ith a more dense, strong and stiff reinforcing 
material in fibrous or whisker form. Hence, high strength-to- 
weight and stiffness-to-weight ratios are readily obtained. In 
electronic applications, circuit boards typically have two or 
more different materials in laminated form. 

During the last decade much effort has been directed 
towards developing and demonstrating the properties of filament- 
reinforced composite materials. Until veiy recently these 
efforts have been oonoentrated on the characterisation of the 
basic materials and the development of information on their 
behaviour. The basio difference between isotropic materials 
and fibre-reinforced composites is the variety of couplings, 
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for example, between inplane loading and out-of -plane deflec- 
tions or between applied bending load and twisting response 
that the latter exhibit. Thus on application of an inplane 
load the laminate deflects transversely, while if the pure 
bending is applied the plate will undergo twisting. The 
couplings between inplane and out-of -plane effects may also 
give rise to warping of the laminate when temperature changes 
take place. The presence of such couplings result in appre- 
ciable reduction in stiffness of the laminate in many situations. 
However, if the plies are arranged symifietrically with respect 
to the plate middle surface, a 'balanced' laminate is achieved 
in v/hlch this coupling does not occur. Although balanced 
laminates are common in current designs, it is evident that 
designs will occur in general in which the laminates are not 
balanced. Por example, exercises in the minimum v/eight designs 
of panels could well shov/ thin plates to be desirable, and if 
these have s fewer number of plies the number of available 
options may be strictly limited and the balanced laminate may 
be an impracticable proposition. Situations v/ill also exist 
where for an efficient structure, a variable skin thickness is 
desired as, for example, in a wing structure when the wing 
thickness is reduced towards the tip. A single ply removed 
from a balanced laminate, may create imbalance and introduce 
the couplings referred to. Consequently, coupling effects are 
to bo eliminated as far as practicable in applications of 
composite laminates. Y/hitney and leissa observed that the 
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degrading effect of coupling for an antisymmetrically laminated 

plate can be minimised by increasing the number of layers in a 

2 

laminate of a given thickness. Bartholomev/ showed that 
selection of suitable ply stacking order enables the elimination 
of coupling and thus enables orthotropy to be achieved in 
bending as v/ell as in in~plane loading conditions. This is 
knov/n as special orthotropy. He assumes each ply to be 
identical in thickness and material properties ai'id the fibre 
directions to lie symmetrically about the or-'thotropy axes for 
tho l£uninate. Having aciiieved inplane orthotropy, the 
remaining requirements for designing a specially orthotropio 
laminate are obtained by controlling the ply stacking sequence. 
It is significant to note tViat tliiclaiess distribution of the 
laminae of a laminate of a given thickness has not been 
exploited to obtain extensional and bending orthotropies of a 
laminate and a part of the present study is devoted to fill 
this gap . 

1 .2 BUCKLIN& OF LAIdlHATED PLATES 

Since the mix-sixties the use of anisotropic fibre- 
reinforced composites as structural members has increased 
considerably. This increase owes to the fact that the designer 
con t^e advantage of the anisotropic properties of the 
materials and thereby use this material quite efficiently. Of 
particular interest are aiiisotropio or layered composite 
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materials in the form of thin plates subjected, to inplane 

compressive loads. Q?he primary mode of failure of these plates 

is buckling. Considerable amount of information is available 

5 4 

on the buckling of structure panels of isotropic materials, ’ 
the buckling of orthotropic panels^ and on infinitely long 
anisotropic plates under special load conditions 

The buckling of finite homogeneous and layered 
anisotropic composite plates has been receiving serious 
attention since 1967. Chamis used the Galerkin method to 

determine the buckling response of simply supported rectangular 
plates under uniform normal and shear loads. The plates were 
assumed to be orthotropic along the material axes which are 
oriented at an angle with the load directions. The buckling 
loads are determined by the power method whioh is applicable 
for the determination of dominant eigenvsiues. Bending and 
stretching coupling effects are also discussed. The complexity 
of the problem is, hov/ever, increased with the application of 
variable load along the edges. To our knowledge, this problem 
has not been discussed to date and the present study is partly 
an attempt in this direction. 

Chiu^^ employed the energy method to determine the 
buckling loads for unstiffened orthotropio plates with a' 
variety of boundaiy conditions and for stiffened plates with 
supported free bounda^ conditions. Buckling mode shapes are 
represented by trigonometric series which, in general, do not 
satisfy all the boundary conditions. The unsatisfied boundary 
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conditions are handled through the use of the Lagrange multi- 
plier method. 

1 5 

Massey considered the elastic buckling of ortho- 
tropic rectangular plates by using the 'Reduction method' and 
obtained the elastic buckling coefficients for flat, rectangular, 
orthotropic plates v/ith uniform inplane loads in one or two 

directions for various types of edge conditions. 

1 6 

Harris presented a detailed discussion of the 
influence of elastically restrained inplane displacement on the 
buckling load of orthotropic rectangular plates subjected to 
end compression. His analysis was based on v/ell-known equations 
for elastic buckling of simply supported orthotropic plates 
subjected to biaxial inplane loading and with no coupling 

between inplane and bending effects. 

17 

Ashton and Waddoups employed the Eitz energy 
formulation to predict the buckling response of anisotropic 
rectangular plates with various boundary conditions. The 
displacements were expanded in terms of beam functions. In 
this , the moment boundary condition was not exactly satisfied 
for simply supported or free edges , but the results were found 
to be acceptable. 

IQ 

Ashton identified the analogy between skewed 
isotropic and anisotropic plates. This allows direct adaption 
of solutions for skewed isotropic plates to rectangular 
anisotropic plates for which and I>2^ are not zero. 
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19 

Wang observed that a separation of variable 

technique cannot be employed to obtain an exact solution for 

simply supported anisotropic plates. The exact satisfaction 

of the moment boundary condition precludes the use of separation 

of variables, if the love-Kirchoff hypothesis is accepted. 

20 

Holston developed an approximate equation for the 

buckling load of orthotropic plates v/ith three edges simply 

supported and the remaining edge free. 

21 

Richer and Mandell have attempted to outline the 

differences between the buckling behaviour of laminated 

composite plates and isotropic plates. They derived the 

governing differential equations utilizing the variational 

approach and including the effect of coupling betv/een exten- 

sional and bending stiffnesses. They have shovwi that a 

reduced flexural stiffness method can be used to account for 

coupling in some types of laminates. However, Reissner and 
22 

Stavsky were the first to identify the full complexity 
associated with the analysis of unbalanced Eingl e-ply 
laminates. They analyzed the response of a two-layered, angle- 
ply composite plate in terms of a stress function and out of 
plane displacement. 

25 

Whitney and Leissa analyzed the antisymmetrically 

laminated angle-ply and cross-ply plates in terras of a stress 

function and out of plane displacement, w. They extended the 

22 

v/ork of Reissner and Stavsky by presenting solutions for the 
transverse vibration and buckling analysis of angle-ply 
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composite plates. The exact solutions for the buckling load of 
an angle-ply laminate with all edges simply supported and 
including the effect of coupling were obtained. However, the 
results for the buckling of cross-ply composites were not 
presented. 

24 1 

IVhitney and VThitney and Leissa formulated the 

analysis of cross-plied and angle-plied composite plates in 
terms of three equilibrium equations in terms of displacement 
variables u, v and w. The stability analysis^ of simply 
supported angle-ply plates is cast in terms of the displacements 
associated with buckling and the loads prior to buckling. The 
buckling analysis of simply supported cross-ply plate is again 
omitted. 

25 

Turvey determined the buckling response of moderately 
thick laminated plates simply supported on all sides. He used 
the concept of reduced plate stiffness along with shear correc- 
tion factor to obtain an approximate solution for biaxial 
buckling of simply supported laminated square plates. Design 
data is generated in the form of buckling load versus thickness 
ratio curves for angle-ply plates fabricated from a high modulus 
lamina. 

Jones'^ considered the buckling and vibration of 
unsyrame trie ally laminated cross-ply rectangular plates simply 
supported on all sides. The classical laminated plate theory 
including the Kirohoff hypothesis of non-deformable normals is 
used for simply supported rectangular plates. 
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27 

Chatnis presented the theoretical results for the 

buckling of anisotropic plates subjected to simple and combined 

inplane loading. The plates are made from fibre composite 

material of boron/aluminium or high modulus graphite/resin. 

The results are presented in non-dimensional form as buckling 

load versus fibre orientation angle for various plate aspect 

ratios. An excellent review of the v/ork done on the mechanics 

28 

of composite materials has been reported by Bert 

A survey of foregoing literature on the buckling 
aspect of composite materials indicates that practically very 
little v/ork has been done to study the effect of the following 
on the buckling response of antisymmetrically laminated angle- 
ply and cross-ply plates: (i) boundary conditions other than 
simply supported (ii) thickness ratio of the laminae of a 
laminate of a given thickness and (iii) cross-ply ratios. The 
purpose of this thesis is partly to fill this gap by attempting 
an efficient formulation to yield closed form solutions to 
predict the buckling response of certain types of laminated 
composites for various boundary conditions and investigating 
the effect of the various couplings and thickness ratios of 
laminae and of oross-ply ratios on the buckling load of 
laminated composites. These results have been put in the form 
of design charts. 
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1.3 BUGICLDTG 01 LAIMITATED PLATES TOPER lATBRAI, LOADS 

It is worthwhile studying the effect of transverse 
load on the buckling of the laminated plates. In effect it is 
like an elastic foundation on v/hich these plates are resting. 

The problem of isotropic beams and plates resting on elastic 
foundations have received extensive treatment in literature?^”"^^ 
However, it appears tha.t onljr a few solutions"^^”^^ have been 
published for the case of anisotropic composite plates on 
elastic foundation. In this thesis, the Winkle r-Pastemak 
foundation has been chosen for investigating the effect of 

lateral loads on buckling of laminates. 

44 

Purvey obtained a simple exact solution for the 
static response of a simply supported, antisymmetrically 
laminated cross- and angle-ply strips resting on WinlcLer- 

Pastemak foundation and subjected to a uniform lateral load. 

45 

The same author in another paper analyzed a uniformly loaded, 
simply supported antisymmetrically laminated cross- and angle - 
ply rectangular plate resting on a general (Winklex'-Pastemak) 
elastic foundation. The importance, of the material coupling 
and foundation stiffness on the plate response has been 
discussed. 

A survey of above literature reveals that a study of 
the influence of Winkler-Pastemak stiffnesses on the buckling 
response of antisymmetrically laminated cross- and angle-ply 
composite plates has not been reported so far. 
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1 .4 SCOPE OF THE THESIS 

The thesis has eight chapters. In the second chapter, 
the characterization of the angle-ply laminates v/ith special 
reference to antisymmetric laminates is considered. It is 
observed that a suitable choice of a thickness ratio of the 
laminae of a laminate can lead to special orthotropy v/ith a 
consequent increase in the overall stiffness of the laminate. 

Chapter 3 is concemed with the formulation and method 
of analysis for the static and buckling response of angle-ply 
and cross-ply plates. A system of three governing equations 
derived from energy a,pproach is reduced to an eighth order 
partial differential equation in terms of a. displacement 
function. The displacements and stress resultants are also 
expressed in terras of a single displacement function. 

In the fourth chapter, a specimen problem entitled 
"static analysis of an antisymmetrically laminated angle-ply 
plate subjected to pure bending with tv/o opposite edges simply 
supported and the other two free" is solved and the results 
thereof in a limiting case are compared with those available in 
literature to confirm the validity of the method proposed in 
Chapter 3« It is shown that the proposed method leads to a 
faster convergence. 

In the fifth chapter, buckling response of antisymme- 
trically laminated aiagle-ply and cross-ply plates with two 
loaded edges simply supported and various boundary conditions 
on the remaining two edges is discussed In detail. 
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The analysis of Chapter 5 is extended in Chapter 6 
to study the influence of V/injk:ler-Pastemak stiffnesses on the 
buckling response of composite plates with various boundary 
conditions. Some of the results .presented herein appear to be 
new and reveal several interesting features. 

Chapter 7 is concerned with the buckling response of 
simply supported antissmunetrically laminated angle-ply plates 
subjected to pure bending at the two opposite edges. The 
problem is solved using Galerkin method and the limiting 
results (with uniform xmiaxial loads) are compared v/ith those 
available in literature. 

The thesis concludes with summing up of the major 
conclusions dra,wn from these studies together with suggestions 
for further work. 


o 
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CHAPTER 2 

CHARACTERISATION OP THE LAMINATED COMPOSITES 

2.1 INTRODUCTION 

In recent years the use of anisotropic laminated 
composites as structural members has increased considerably. 

This is due to the fact that these can be used yery efficiently 
by taking advantage of their anisotropic properties and light 
weight coupled v/ith high strength. A considerable amount of 
research work, therefore, is under progress regarding the 
elastic behaviour of laminated composites. Until very recently 
these efforts have been concentrated on the characterisation of 
the basic material properties and the documentation of infor- 
mation on their behaviour in practical applications. 

The laminated plates consist of a number of layers 
of thin ortho tropic sheets v/ith the ortho tropic axes of symmetry 
of each lamina oriented at an arbitrary angle to the plate axes. 
Two of the important types of arrangements are (i) antisymme- 
trical cross-ply and (ii) antisymmetric al angle-ply laminates. 
Antisymmetrical cross-ply laminates consist of an even number 
of layers of equal thickness and same elastic properties with 
the orthotropic axes of symmetry in each layer alternately 
oriented at angles of 0° and 90° to the plate axes. Antisymme- 
trical angle-ply laminates also consist 6t an even number of 



13 


layers having the same thickness and elastic properties but the 
orthotropic axes in each ply (layer) are alternately oriented 
at angles at + © and - 0 to the plate axes. It is v/ell -known 
that antisymmetrical cross-ply and angle-ply laminates exhibit 
a variety of couplings, for example, betv/een inplane loading 
and out-of-plane deflection or between applied bending load 
and twisting response. The presence of such couplings in 
addition to complicating the analysis of laminated plates, 
result in an appreciable reduction in overall stiffness of the 
plate for certain types of laminates. Consequently, it is 
desirable that' coupling effects are eliminated as far as 
possible in the application of composite laminates. 

A laminated plate so designed that the coupling 
between inplane and bending effects does not occur and which 
is orthotropic both inplane and in bending is generally 
referred to as ’specially orthotropic ' . The analysis required 
for such plates is more straightforward than for more general 
laminated plates. 

A tevf methods have been suggested by various investi- 
gators to do away with the degrading effects of couplings, In 
this chapter, the various methods for eliminating the adverse 
effects of couplings associated with the general laminated 
plates are reviewed and a technique is suggested to arrive at 
the decoupling effect-. It is further shown that in the case 
of angle-ply laminates for a ttacking sequence of (45° + 0), 
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the thickness distribution of laminae at which the decoupling 
effect is achieved, is independent of material properties. 

2 . 2 KEVIEW OF LITBRATUBE 

The effect of coupling between inplane stretching (or 

membrane action) and plate bending as discussed above was 

46 

recognized as early as 1953 by Ambartsumyan who was concerned 

with the study of unsymme trie ally laminated orthotropic shells. 

47 

At about the same time, Reissner also realized the presence 

of coupling in laminated plates. This coupling effect is the 

major difference between the macroscopic structural behaviour 

of an arbitrarily laminated plate (either isotropic, specially 

orthotropic or anisotropic) and homogeneous plates (either 

isotropic, specially orthotropic or anisotropic) such as 

48 

treated in books by Timoshenko and Woinowslcy-Krieger and 
6 

Lekhnitskii°. 

49 

Werren and Uorris have sho\m that it is possible to 
orient the layers of inultiple-layer laminates in such a way 
that the resulting elastic stretching behaviour is isotropic 
i.e., it has elastic coefficients which are independent of 
laj'-er orientations in the plane. According to them, the 
conditions which must be met by the Imninating sequence to 
achieve this result are (i) the total number of layers must be 
three or more, (ii) the individual layers (denoted by index k, 
ranging from one to n) must have identical orthotropic elastic 
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coefficien’ts and thicknesses and (iii) a typical layer k, must 

be oriented at an angle = it (k - 1)/n. with respect to a 

reference direction. However, the laminate made according to 

40 

the Werren-Norris design is isotropic in regard to stretching 
and not, in general, in regard to bending and stretching -bending 
coupling, 

50 

Abrams and Scheyhing presented numerous configur- 
ations for which the elastic behaviour of the laminate is 
identical to that of a homogeneous orthotropic or isotropic 
material. Herein, the constitutive equation of a thin multi- 
layer plate or shell element is derived and compared with its 
counterparts for homogeneous media in order to determine the 
mathematical relations which will ensure composite orthotropic 
and isotropic elastic behaviour, laminate configurations 
satisfying these relations are then found by approximate and 
exact means, 

2 

Bartholomew has shown that selection of a suitable 
ply stacking order for uncoupled, orthotropic laminates enables 
orthotropy in bending, as well as in-plane, to be achieved, 

This paper assumes some degree of standardisation of design 
parameters. Each ply is taken to be identical in thickness 
and material properties, and the fibre directions lie 
symmetrically about the orthotropic axes for the laminate. 

Having achieved in-plane orthotropy, the remaining requirements 
for designing a specially orthotropic laminate are achieved 
by controlling ply stacking sequence. 
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From the above it can. he seen that not all the 
parameters have been considered to achieve uncoupled ortho tropic 
laminates, An attempt is, therefore, made in this thesis to 
exploit the thickness variation of laminae constituting a 
laminate of a given overall thicloiess to achieve complete 
elimination of coupling effects (special orthotropy) for anti- 
symmetrically laminated angle-ply plates. An analysis leading 
to partial reduction of coupling effects for more general types 
of laminates is also made. It has been further shown that for 
a stacking sequence given by (45° + 0), the thickness ratios 
of the laminae at which the complete elimination of coupling 
i terras occurs are independent of the material properties. 


2.3 ANALYSIS 


Consider a laminate with N layers (Figure 2,1), Each 
ply of the laminate, made of fibre reinforced material, is 
assumed to be orthotropic with uniform material properties E.^ , 
Eg, G ^2 ®*3.d where 1 and 2 refer to the axes along and 

perpendicular to the fibre. Following the usual notation as in 
reference 2, the matrix of elastic stiffness coefficients for 
such a ply is given by 




Q|2 

0 


Q 


12 


Q 


22 

0 


0 1 

0 

^66 


( 2 . 1 ) 
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where 


E 


. '=11 


1 


Q 


Q 


'(1 - ^ 2 ^ 21 ^ 




■^21^1 


12 


22 


(1 ~ ■^2l) 

- ^2 

( “ '*^12 "^ 21 ^ 


rr^ 


12 ^'21 


T 


^66 ~ ^12 


The Poisson's ratio -02^ can be obtained from the reciprocal 
relation 


^12 


'21 


(2.2) 


The matrix of elastic stiffness coefficients (2.1) is referred 
to the orthotropic axes. However, if the fibres of such a ply 
are taken to lie at an angle 0 relative to the x-axis of the 
reference plane (Figure 2.1), then the transformed ply stiffness 


is 

given by 





Si 

S 2 

Se 

Q 


S 2 

S 2 

Sg 



..^16 

Sg 

%6 


(2,3) 


where 
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cos^Q + 2(Q^2 + 2^66^ Bxn^9 cos^O + Q 22 sin^O 

^12 “ (Q-ll + ^22 *" ^^66^ sin^O cos^O + Q^ 2(^11^^® + cos^O) 

Q22 = Q-jl sin^O + 2(Q^2 '*■ ^*^66^ sin^0 cos^O + Q22 cos"^© 

^16 ” (Q-i^ ~ Qi 2 - 2Q65) sin© cos^O + (Q^2 ” *^22 

2Qgg) sin 0 cos© 

Q25 = (Q-ii - Q-12 " ^Qgg) sin^© cos© + (Q^2 “ “^22 

2Qgg) sin© 003^6 

^66 ^ (^11+ ^22 *" ^*^12 ~ ^^66^ sin^© cos^G + Qgg(sin^0 + 

cos"^©) 

(2.4) 

A ply with fibres at an angle -. © to the x-axis may 
be seen to have the same elastic stiffness coefficients as 
given in relation (2,4) above apart from Q^g and Qgg which 
change sign. 

Uow, assuming small deformations and following the 

thin plate theory, it can be shown that the relationship 

characterising the response of laminates to applied membrane 

51 

stresses and bending moments is linear, A full 6x6 

symmetric matrix is reqjuired, in general, to describe this 

51 

elastic behaviour. It is given by 


fN 1 

Ta 

\ 

Bl 

i 

1 s^-r 





^ \ \ 

(2.5) 

M 



. [^) 
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where 

N = (H^, Ny, If^) 

M = My, M^) 

are membrane and bending stress resultants and 

£0 = (eo.eo.eo ) 

' X ’ y’ xy' 

K ~ '"^y* '^'xy^ 

are middle plane strains and curvatures respectively. The 
matrices A, B and D are defined as follows; Assume the kth 
ply of the laminate as shown in Figure 2.1, to lie between 

and Zy. Then the components of the membrane stiffness 

matrix A, coupling matrix B and bending stiffness matrix D ore 


given by 

a sinomation 

over the IT plies; 



H 



li 

(4 - 4-1 > 

(2,6) 


= 1 f 

*4 ■ 4-1^ 



It can be shown that A,|g = A^^ 

= 0 leads to in- plane 


orthotropy; B-jg = Dgg = 0 ensures orthotropy in bending while 
^ = 0 yields; an uncoupled lay-up . 

In v/hat follows, the symmetry of the laminate about 


its middle plane is important since the plies considered herein 
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are grouped in pairs, symmetrically placed about the roid-plane 
and numbered outv/ards from the centre (Figure 2.1). When both 
plies of a pair have fibres aligned in the same direction, the 
pair is referred to as symmetric, the term antisymmetric is 
applied to pairs with fibre directions of + 0 and ~ 0 respec- 
tively . , 

In general, it has been noticed that antisymmetric 
laminates tend to become orthotropic for N >- 8. It should, 
therefore, be our endeavour now to investigate v/hether it is 
possible to obtain orthotropic and uncoupled behaviour even in 
cases where W < 8 by manipulating the thickness distributions 
of layers. The minimum number of layers in such a case will 
have to be four and accordingly the analysis is given for 4 
and 6 layers in the following sections. 

2,3.1 Four-Layered Angle-Ply Laminate 

Consider a four-layered laminate with laminae thick- 
nesses t^ and ^2 on each side of the middle siirface and v/ith 
fibre directions aligned at + 0^ and - 02 to the x-axis of the 
coordinate system as shown in Figure 2,2, 

The first of Eq, (2,6) for the present case is 
reduced to 

hj = 2Qi3(ep tj] (2.7) 

where 
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a . . 

13 




aJid Q^^(©) are the reduced stiffnesses evaluated at angle of 

1 J 

orientation of the fibre, 6. 

Por an antiayiranetric angle-ply plate , a . . = 1 and the 

X J 

extensional stiffnesses A. . reduce to 


^ij " ^ij 

for 0/^6 

(2.8a) 

= 0 

for 3=6 

(2.8h) 


where t is the overall thickness of the laminate* 

The relation ( 2 *83) shows that the extensional 
stiffnesses A. . are independent of the individual ply thiok- 
nesses. 

The second of Eg. ( 2 * 6 ) upon substitution of parameters 
assumed in Figures 2,1 and 2.2 and after simplification reduces 
to 

= I ^ ^3 - 4 '> * 

= [st| - 2{t^ tg)® + 2«„t2] 

= + “13) - 

or 

f(t^, t^, a^.) 


- 2z 2 + z|)] 


( 2 , 9 ) 
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The purpose of the present investigation is to reduce 
B . . to zero in order to obtain uncoupled lay-ups. B . • can be 

' Xj 

zero, provided f(t^, t 2 > zero since cannot be 

equal to zero. This leads to 

+ 2R^ - = 0 (2.10) 

where 

R^ = t^/t2 

Eq. (2,10) is quadratic in R^ and yields 

R^ =. - 1 j- fTTo^ (2.11) 

Por an antisymmetric angle-ply laminate, = 1. Hence 

R^ = 0*4142 (2»12) 

The bending stiffnesses from the third of Eq. (2.6) can be 
reduced to 

3 

^ij " for D 6 (2.13a) 

= 0 for 3 = 6 (2.13b) 

which are obviously independent of the thickness of an 
individual lamina for the antisymmetric angle-ply laminate. 

2,3.2 Six -Layered .Angle-Rly laminate 

Considei^ a six layered angle-ply laminate with thick- 
ness distribution and laminae orientation as shown in Figure 2.3* 
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The extensional stiffnesses A. . and bending stiffnesses 
Dj . can be derived as before and are given by 

X J 


hi 

= 2Q^.(0^) [t^ + 

] for 2 ^ ^ 

(2.14a) 


= 0 

for 3=6 

( 2 . 14 b) 

hi 

3 

= 1 Qij(ei) [|- - (•^2 + ■^3)^ + 

“i3 \- *3 (^2 + 



"Pij *31 ■ 

for 3 6 

(2.15a) 


= 0 

for 3=6 

( 2 . 15 b) 

where 




Pij = Qij(e3)/^Ti.(e,) 



For 

an antisymmetric angle-ply laminate, 



“ij Pij ~ 



and 

that 

upon substituting these in Eqs, 

(2.14) and (2.15), 

we find 



t^) 




for 3 7^ 6 

( 2 , 16 a) 


= 0 

for 3 = 6 

( 2 . 16 b) 


^ij = f ^ij^®i^ 





for 3 6 

( 2 . 17 a) 


= 0 

for 3=6 

(2.17b) 
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Eqs . (2.16) and (2.17) shov/ that A . . and D . . are 

J J 

independent of the individual ply thicknesses. 

The coupling stiffnesses for the six-layered 
configuration as obtained from the second of Eq* (2*6) are 
given by 




= Qj ^^(©- j ) {.{^2 ^ 3 ^^ "* ^^1 ^2 + 

+ {(±2 + t,)® t2] 

= ai3(9i) f(ti, tg, tj. 


(2.18) 


B,. = 0 


if f(t^, t^, Pj_^) = 


0 


This gives upon algebraic simplifications 


aj_. rI + 2{a^.- R^)R2 “ (R^ + 2R^ + Pi-j) = 0 


(2.19) 


where 


= t-j/t^ 
^2 “ 


In teims of , Rg can be written as 


E, 


“iD 


( 2 . 20 ) 
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I'or an antisymmetric laminate, Eq^, (2.20) reduces to 

R 2 ^ (R^ ~ 1) + + 2 (2,21) 

2.4 RESULTS MR DISCUSSIOH 

The variations of coupling stiffnesses and B 2 g as 
a function of thickness ratio, R^ (= t^/tg) of the laminae of 
a four-layered antisymmetric angle-ply laminate are shown in 
Figure 2.4 for three different fibre orientations 50°, 45° and 
60°. It is seen that in each case B^g and B 2 g vanish simul- 
taneously at R^ = 0 . 4142 . For other values of R^ , there exist 
finite values of B^g and Bgg. The coupling stiffnesses do not 
exist for extreme fibre orientations of 0° and 90° for obvious 
reasons. Thus it is possible to eliminate the coupling terms 
entirely and achieve special orthotropy at any fibre orientatior 
of a four-layered antisymmetric angle-ply laminate. However, 
if the laminate is not antisymmetric i.e., if the angle of 
orientation of the fibres in each layer is different (+ , 

- ©2 5 + “ ®1 example), it is possible to 

eliminate one of the coupling stiffnesses B^g and B 2 g. Such a 
choice may be used for the elimination of more predominant 
coupling stiffness. Table 2.1 gives the values of a^g and 0:25 
for various stacking sequences for laminates made of different 
materials. These values can be used in Eq. (2.11) to obtain 
a thickness ratio which results In the elimination of one of 
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the coupling stiffnesses orB^g). Figures 2,5 and 2,6 

shov/ the variations of the coupling stiffnesses and B^g 

with Rg (= for given values of (= t^/t^) of 1,0 and 

2,0 respectively for a six-layered antisymmetric angle-ply 

laminate. Here again, a unique value of thickness ratio exists 

for which the coupling terms vanish at any fibre orientation. 

Figure 2.7 shows the plot of R^ versus R 2 for a six -layered 

antisymmetric laminate. This is valid for any fibre orientation. 

A choice of point falling on this curve ensures the absence of 

coupling terms and consequent achievement of special orthotropy. 

The intersection of this curve v/ith the Rg-axis (R^ = 0) yields 

the critical thickness ratio of 0.4142 for a four-layered 

laminate, since with R^ =0, ’the six-layered laminate is 

analogous to a four-layered one. However the case of t^ = 0 

which also corresponds to one of four-layered cannot come out 

as a limiting case of the above derivation since all the 

quantities tend to infinity at different rates. 

Similar to the case of four-layered laminate, the 

absence of both the coupling stiffnesses mey not, in general, be 

< 

possible if the fibre orientation in each lamina is different. 

In such a situation, Table 2,1 may still be used for a particular 
stacking sequence to obtain a partial reduction in coupling. 
Consider, for example, a six-ply laminate with stacking sequence; 
30o/-45‘'/60o/”60o/45°/-30° . 
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The values of a. and P • • to be used in Eq. (2,20) as 
obtained from Table 2.1 are: 

= 0.82385 , = 0.42696 

«26 " 1 ‘92957 , ^26 2.34212 

A choice of and R 2 leading to partial elimination 
of coupling stiffnesses (B.jg = 0 or Bgg = O) for a six-layered 
laminate with the above stacking sequence is given in Figure 
2,8 for a carbon-fibre reinforced plastic (CPRP) laminate. The 
intersection of the two curves which ensures simultaneous 
elimination of and Bgg is seen to be a peculiar behaviour 
of this stacking sequence. It is further observed that the 
point of intersection namely R^ =0,25 approximately and 
Rg = 0,60 approximately is independent of material properties 
as shown by Figures 2,8 to 2.10 for three different materials. 

A similar behaviour, in general, can be expected for a stacking 
sequence of ('n;/4 + 6)» This peculiar behaviour may not, in 
general, be expected from an arbitrary stacking sequence. 

Figures 2.11 and 2,12 drawn for two different materials (G-FRP 
and CPEP) between and Rg for a stacking sequence of 45^/~30o/ 
60°/-60°/30°/-45° , for example, do not exhibit any intersection. 


2.5 COHOIIJSIOITS 

1. For antisymmetric angle-ply laminates, it is possible to 
eliminate the coupling effects by a suitable choice of 
thicknesses of laminae. 
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2 , The coupling effects can be alleviated for general angle- 
ply laminates by a suitable choice of thickness ratios. 

3, It is possible to obtain the orthotropic behaviour of an 
angle-ply laminate with a fewer number of layers 

(4 ^8) by suitably selecting the laminae thicknesses. 

- o - 
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TABLE 2.1, VALUES OE oc^ g AND VARIOUS STACKBTG SEQUEUGES 

EOR CORIMOU COICPOSITE MEEIABS 

I. Carbon Fibre Reinforced Plastic 

(E^/E^ = 6.333j G^2/®2 = 0*667} ->>^2 == 


Sr, Me 

1 

t 

. } Stacking Sequence 

! - 

1 1 

' oc ' 

! 16 } 

_ ! _ - - 1 

“26 

1 

30oA45°/45°/“30o 

0.8238 

1.9295 

2 

45V-30o/30oA45° 

1.2138 

0.5182 

3 

30°/-60°/60o/-30o 

0.4269 

2,3421 

4 

60o/~30o/30o/-60o 

2.3421 

0.4269 

5 

60V-45°/45o/~60^ 

1.9295 

0,8238 

6 

45°/-60o/60o/-.45o 

0.5182 

1.2138 

II. 

Graphite Fibre Reinforced 

Plastics 



(E/E2 = 40.0} G^2/®2 

0.50; ^^^2 = 0.25) 


1 

30%45°/45V-'30® 

0.7731 

2.2805 

2 

45%30o/30o/~45o 

1.2935 

0.4385 

3 

30o/~60o/60o/~30o 

0.3390 

2.9496 

4 

60o/-30o/30'>/-60o 

2.9496 

0,5390 

5 

60 °/^ 45 ° A 5 °/- 60 ° 

2.2803 

0.7731 

6 

45V’-60V60o/-45° 

0.4385 

1.2937 


Continued, . * 
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TABLE 2,1 (continued) 

III. Boron Fibre Reinforced Plastic 


(E/Eg 

= 10.0; ^ ^>.333? 

A 2 = 0-30) 


" 1 

Sr. No. { 

r 

1 

f 

Stacking Sequence j 

- - t ... 

' r 

«16 I 

- - - -1 

oc^ ^ 

26 

1 

30V-45V45V-'30o 

0.7676 

2.3298 

2 

45V-30°/30o/-45® 

1.3028 

0.4292 

3 

30oA60o/60o/~30o 

0.3294 

3.0353 

4 

60o/-30°/30o/-60o 

3.0353 

0.3294 

5 

60oA45V45°A60o 

2.3298 

0.7676 

6 

45°A60o/60oA45° 

0.4292 

1.3028 


IV. 

Glass 

(E/E, 

Fibre Reinforced Plastic 

•p — 3.0 5 g^^/e^ — 0,505 

= 0 . 25 ) 


1 


30oA45°/45°A30o 

0.8334 

1.8788 

2 


45o/„3oo/3ooA45® 

1.9980 

0.5322 

3 


3 O 0 A 6 O 0 / 6 O 0 A 300 

0.4436 

2.2542 

4 


60oA30°/30oA60« 

2,2542 

0.4436 

5 


60°A45°A5oA6O'> 

1.8788 

0.8334 

6 


45V-60V60oA45° 

0.5322 

I. 1998 







FIG. 2-E: GEOMETRY OF A FOUR LAYERED 
LAMI NAT E. 





FIG 2-4 VARIATION OF COUPLING STIFFNESSES W(TH 
THICKNESS RATIO FOR A FOUR-LAY £ R £ D 
A N T I S Y M METRIC LAMINATE 
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FiG 2 5 VARIATION OF COUPLING STIFFNESSES 
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FIG 2 7 LOCUS OF ORTHO- 
TROPIC THICKNESS RATIOS 





ELIMINATION OF COUPLING STIFFNESS 
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FIG. 2 9 LOCUS OF THICKNESS RATIOS FOR THE PARTIAL 
ELIMINATION OF COUPLING STIFFNESS 
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FIG2ti LOCUS OF THICKNESS RATIOS FOR THE 
PARTIAL ELIMINATION OF COUPLING STIFFNESS 
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CHAPTER 3 

EORI^ULATIOH OP TliE PROBLEM AND L1ETHOD OP SOLUTION 


'3 • 1 I NTRODUCTI ON 

Laminated plates are one of the si 2 nplest and most 
vvidRly used practical applications of composite laminates. 

There exist many kinds of theories and levels of sophistication 
in the theory of laminated plates, both liiiear and nonlinear. 
These can be subdivided ass (1) thin plate analysis (2) moder- 
ately thick plate analysis including the thickness-shear and 
thickness-normal deformation (3) cl-assical elasticity analysis 
applicable to non-homogeneous laminated plates. Ov/ing to 
simplicity of analysis, thin plate theory has been most popular. 
The simplified form of this theory has been applied to a plate 
laminated symmetrically v/ith respect to its middle plane. ^ 

A general form of this theory has been applied to arbitrarily 
laminated plates including bending-stretching coupling 
effects.^^’ A more complete Glassification of laminated 
plate problems is found in reference 55. 

The response of laminated plates is governed by three 
coupled displacement equations of equilibrium. The laminated 
plate analyses have been carried out mostly by assuraing 
suitable functions for the three displacement components 
separately. However, to obtain a more general solution, it is 
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dosirable "to combine the three displacement ec^uationa of 
equilibrium into a single higher order partial differential 
equation, Ih what follows) a system of three coupled goyeming 
equations for arbitrarily laminated plates is derived by an 
energy approach and this system of three equations is then 
reduced to a single eighth order partial differential equation 
in terms of a single displacement function. A closed form 
solution to this eq,uation is then found for some types of 
laminates under various boundary conditions, 

3.2 DERIVATION OE EQUILIBRIUM EQUATIONS AND NATURAL BOUNDARY 
OTpITTQlf^ ' 

3,2.1 Assumptions 

( 1 ) Each layer is orthotropic , The principal material 
directions of each layer are, in general, aligned at an 
angle 9 to the plate axes. 

(2) The plate thicioiess is very small compared to its length 
^d width (plate is assumed to be thin) . 

(3) The usual assiMptions of the small deformation theory of 
thin plates are valid here also, 

(4) Strains e^, Sy and are small, 

(5) No body forces exist. 
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3.2.2 . Strain-Displaoement^Gurvature Relationships 

The rectangular coordinate system considered is shown 
in Figure 3.1. The assumptions of classical thin plate theory 
lead to the following displacement relations: 

u = u°(x, y) - z || 

V = v°(x, y) - & 0 (3.1a) 

u 

w = w (x, y) 


The corresponding strain-curvature and the strain-displacement 
relations are 




e° + z ir^ 




®x ~ 







(3.1h) 


z 




where 





z 

0 

- ^ , 

?JX ' 

0 

-0 - . 
y ^7 ' 

0 

^xy 

dj- -3x 




’ ^'xy 

-2 

" " ’ 

= - 3-7 

fix <^y 

The superscript ' o • 

refers to the 

middle plane 

values. 


3.2.3 Derivation of Equilibrium Equations and Boundary 
Conditions by Energ^ Approach 

A general expression for the elemental strain energy 
is given by 
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1 

2 


dV = 7T + cv..e_ + T 


XX y y z 2 
dx dy dz + q(x, y) w dx dy 


e +T e +i; e ], 
2cy xy yz ys xz xz 


(3.2a) 


Incorporating the general assumptions of classical thin plate 
theory and integrating over the thickness t, the strain energy 
of the plate is given by 


_ iff h] 'h^ + il- + p 


^ ~ 2 i \ ‘'"x d X 

"" Ik- 


y ay 'xy ^dy dx 

'•2 2 2 

(f w <7 


T,/r r V/ TVT (f W r,n/r <7 W n n , 

“x ^ \ ^ - 2^ .Tay ay 

+ f ( '3(x, y) w ax dy 


where 




X 


t/2 

■ j O' dZ5 
-t/2 ^ 


(3.2b) 


t/2 t/2 

"y = -J/2 V*>- "xy = ://, 


xy 


t/2 

= j z dZ 5 M, 

-t 72 


V2 t/2 ^ , 

y = j' ^y 2; dZ5 M = r '^xy 

^ -t/2 y -t/2 


are the inplane normal and shear forces (N , U , l'T„„) and 

^ y ^y 

bending and twisting moments (M , M , M„ ) respectively. 

tj "^y 

55 

The constitutive relation is given by 


'"x 

1 

Ai ^ 

^12 

■^16 I 

®11 

®12 

®16" 


r 1 

X 

"y 


; -^12 

^22 

•^26 « 

I 

®12 

®22 

®26 


e° 

y 

"xy 


■^16 

^26 

•^66 ‘ 

^16 

®26 

^66 


e® 

xy , 

Kx 

I 


B 11 

h2 

®i6 ; 

^11 

®12 

“16 


I 

i 

•v 


®12 

®22 

®26 ’ 

I 

D-,2 

^22 

^26 


ICy : 


t 

i 

i 

j 

_®16 

®26 

^66 ;; 

^16 

^26 

^66 

I 

i 

1 

^xy J 
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v/here 


A. . , B . . , D . . = 

^10 » -“ij’ 


t/2 


:(k) 


r Q>7 (1, 2, z-") dz 
-t/2 ^ 


(3. 2d) 


““ f k ') 

Q> is the reduced stiffness of the kth lamna. fhe 
substitution of relation (3.2c) into Eq, (3.2b) yields 


1 i i 

2 i 


'•up 2 o n o o 2 


3 O T O 

+ PA iiL. . ^ 

■^16^ 3x ?.y 


p 0 'j o n 0 " O O O 

pu 5 . /lv_ V , ^U dv N 

^ Y ^ ■^26^ Sx 3y oy 3y ^ 




u 


66'‘lTy~' 


^ o 2 - 

ai"^ “ IPc 


^2 '> o 

o V/ ou 


■>2 o 0 

1 V/ P V 


2 , 


12/^ 2 ay 


^ 2 


Pu \ 


pp p-p <r9u° /? d^p 

2^22 ^ ^ ^ 


- 2B ^ w p V . “"w C' u 

2 .-,0 


' 


oTi fCLJL ^ 


■) - 4B 


'd U' 

2 X 9 y ?> X 

a 


■) 


66 fix 3 y 


•M 

•9y 


* 21-12 f^T jT? + 


0 

/; Y 


i2 

0 w ^ 


w 


- 2 2 
+ D_(f^) 


- «16 I? ieV ^ 

(3.2e) 


'y 


Considering the variation of each term in Bq, (3,2e) and 
applying the Principle of Virtual Displacement, the following 
s et of equilibrium equations and the associated boundary 
conditions are obtained; 
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Equilibrium Equations 


^>2 o „ 

A p u , A 'V 

^11 ^^2 ^12 ax ay 


>',2 o q2 0 ',2 o q2 o 

c' V ,4 ^ . SLjLa j. a ° V 

axay ^xTy ^^ 2 ^ ^26 3 ^2 


/ a^u° , ^v° \ p 3 p 

1^^ “ ®11 ^ “ ®12 


<3 \{ 

- 5Bi6 


„ B ^ 

rV 


2B 


3 ^w 


‘66 T n 2 

<'.'x c'y 

2, o 


= 0 


(3.2f) 


2 0 i.2 o /j2, o '^2 q q2, o 

^12 axay ''^22 ;jy2 ^^16 g^2 ^26^^ bx Sy + ^^2 

+ + W - ®i2 ®22 3y3 ®16 3x5 


•) 


5B . " .. .w ,, 

^ 3x 3y2 


^ ^w pB - •■ 

~ “ ^-^66 a 2^) 
f/x y 


35 ,, 


(3.2i 


D 


j- 9T1 d ^ W p 3^VV .p 3 ^W 

'11 E? " " “22 e7 " ®16 s,J3y 


2 " ®11 "3 

r d X 


0 3 o 
o u 


_av_) _ B 2ii - 3B 3 5u° 

dx^dy 3y^ 3x^3y 


: B r-i^ + 

3x?y^ dx^ dy 


3 ^v° p / , 


3 V 


A 3 o 


- Bi 6 --g^ --112615 g^gy. „y 


o 3' 

O ■y-' 


? 30 o 3 , o 

- 2Bg6 - 2Bgg -*-^ 

° 3x 3y oxd y 


q(x, y) 
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Boundary Conditions 

Different types of simply supported and olatnped edge 
boundary conditions are possible for laminated plates • The 
detailed discussion of these boundary conditions is given in 
references 56 and 57. The follov/ing types of edge boundary 
conditions (Figure 3.1) are considered in the present thesis s 

(a) Simply supported boundary conditions: 

(i) Antisymmetric angle-ply laminates 

w = 0 , = 0 } 

X = 0, a (3.2i) 

U° = 0, = 0 ,) 

(ii) Antisymmetric cross-ply laminates 

w = 0., \ = 0 J) 

/ X = 0, a ( 5 * 23 ) 

v° = 0 , ^ ^ 



(3.2k) 


( 3 . 21 ) 
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, - B H ' B 

■31 - -®11 ?,3 -®12 


3B 


P;3 


~ B. 

^X“oy f)y' 


1 6 n 2 '■* ~ 26 3 


- 2B 


35 


66 


3 x Sy"" 


132 


a 2 93 

“ B-t o - I ^ ^ Boo ' " ^ ^ B 


95 


- 3B, 


35 


12 a 2'i 22 -A 3 ~ 16 ''■' 3 ” •^’^26 ''• '^■" 2 

ax dy dy^ dx^ c’xdy 


2B 


;)3 


66 o 2 *•■, 

o-x dy 


J 33 


a"*- 

;^ + 2D 
dx 




a ^ 

+ B^o r r + 4D. 




11 B..4 --12 .^2 3y2 -22 3^4 ^^^16 


9 4 

■*■ ^^^26 -i a' 3 
dx dy-" 


3^ 


ax^ay^ 


Assume 

u° = - l ^2 y) “ ^13 y) 

v° = B (x, y) 

w = 1-11 'Mx, y) 


(3.6) 

(3.7) 

(3.8) 


where the functions ti(x, y) and 7'’(x, y) are continuous in the 
domain of the plate. These representations are chosen to 
satisfy Eq. (3*3) identically. 

Substitution of Eqs. (3. 6-3.8) into Eq. (3.4) leads to 

(L22B-|-j ^"12^21^ B (xj y) + (L22B.11 — ^13^21^ (xj y) = 0 


This can be rewritten as 


PB + Qy 


G 


(3.9) 
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If we now choose 

r) = - Q ^ (x, y) 

y = P cp (x, y) 

Eg.- (3.9) is identically satisfied. 

This enables us to express relations (3. 6-3,8) in 
terms of as 


u° 

“ (lj'j2^23 


y) 

( 3 . 10 ) 

v° 

= ( 1131*21 


y) 

(3.11) 

w 

= (L 11 I 22 


y) 

( 3 . 12 ) 


where 

0(Xf y) = ^ (x, y) 

substitution of expressions (3.10-3.12) into Eq. (3.5) leads to 

^^ 11 ^ 22^33 “ ^ 21 ^ 12^33 “ ^ 32 ^ 23^11 ^ 31 ^ 12^23 ''' ^ 32 ^ 13^21 

- ^ (3.13) 

Expanding the differential operators, the above equation can be 
written as 
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where 

Al = ® 11 ^ 2 A^ 6®16 - AggB^^) - + I)^>,(A^^Agg + A^g) 

A 2 = 2 B^^(A^gG + E^B^g - 3 AggB^g - A 2 gB^^) 

+ 2 B^g( 2 A^gB^g - A^^G) + 2 D^^(A^gAgg + A^^A^g 
“ A-|gE^) + 4I)^g(A^^Agg - A^g) 

A3 = ^^II^^A^gBgg + E^G - 5 A 2 gB^g - AggG) 

+ 2B^g(2A^gG + 3E^B^g - 5AggB^g - 3A^^B2g) 

+ D^^(2A^gA2g + Agg + A^^A22 -' ) 

+ 8D^g(A^gAgg + A^^A2g - A^gh^) - ^22^11 "" A-^-^G’ 

+ 2P(A^^Agg - A^g) 

A4 = 2B^^(A^gB22 + 3 E^B 2 g - AggBgg - 3 A 22 B^g - A25G) 

- 2 B 22 B-j 5 A-j 1 + 4 B^g( 2 A^gB 2 g + 2 E^G - 3 A 2 gB^g - 2 AggG) 

“ ^®11^A22At6 ^ ^ 26^66 " A 26 ® 1 ^ 

+ 4I>^5(Agg + A^^A22 + 2A^6A26 “ ''■ ^^26^A-j -jAgg 

~ A^g) + 4 E(A<|gAgg + A^^Agg - A^gE^) 

Ag = 2B^^(E^B22 + A 25 B 26 ~ A22^^ ^^22^Ai5^-lS “ ■^i-jG") 

- 4 Bi 6 A 26 « - 4 B 26 Ai 6 f^ - 9 (A 22 B 2 g + A^^B^g) 

+ 4 B^gB 25 ( 5 E^ - ^Agg) +E^^(A22Agg -- A^g) 

+ B22^A^lAgg - A^g) + 8D^g(A22A^g + AggAgg - AggE^ ) 

+ SDggCA^gAgg + A^^Agg - A^gE^) + 2G(E^ - AggG) 

+ 2 E(Agg + A^^Agg + SA^gAgg - ) 
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Ae ^-^26^22 " ’^22^26^ 

+ 2B22(3E^B^g - A^gG - 3A^ ^B2g - AggB^g) 

+ 2B^g(4A2gB2g -- 3 A 22 G) + 4B2g(2E^G - 2AggG - 3A^gB2g) 
+ 2D22(A^gAgg + A^^Agg - A^gE^) 

+ 4I)^g(A22Agg “ A2g) + 4D2g(A^^A22 + ^A^gAgg + 2Agg 
“ 2®-)) + 4E(A22A-|g + AggAgg ~ A2gE^) 

^7 = ®22^*^-^"26®16 " ^^^16^26 “ " -^11^22^ 

" ^^16®26'^'22 ■'■ " ^^66^ 

+ ^22 1^22 ■'■ -^"66 ~ 

+ 8D25(A22A^g + AggAgg - A2gE^ ) + 2P(A22Agg - Agg) 

+ G(4A2gB2g - A22Cr) 

Ag = 2B22(A2gG + B2gE^ - ^^66® 26 “ ^^16^22^ 

+ 2B25(2B2gA2g - GAgg) + ^^22^^22'^16 '^26^66 ” -^26^1^ 

•i- 4I>26(^‘-66-^22 " -^26^ 

Ag = ®22^ ^-^26^26 *■ ®22-^66^ “ -^22^26 

+ B22(A22A55 >• Agg) 

= A-|2 Aggj E = D^g + 2Dgg| 6 = B^g + 2Bgg 

Equation (3,14) is the governing equation for laminated plates, 
j It can be easily verified that this equation reduces to the 
classical thin plate equation for an isotropic case. Jn v/hat 
follows, we discuss two special cases. 
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3.3.1 iuitisyimnetrlc Angle-Plv Plates 


It has been shown in Chapter 2 that for an anti- 
symmetric angle-pLy laminate 

■^16 ^ -^26 ^ ^ ®12 “ ^22 ~ ®66 ®16 " ^26 “ ° ( 3 . 15 ) 

For such a laminate, the coefficients of the governing Sq, 
(3,14) assume the following values: 




A 


A, 


4 “ -^6 “ ^ “ 

2 

^11 ^""66^11 “ -^16' 

.2 


(3.16) 

(3.17) 


A. 


^11^-^66^11 “ ®16^ 

2B2g(3E^ - 5 Agg) H- 2A^^(AggF - 33^^335) 

+ D^^(A^-|A22 + A55 E^) 

"^^1 6^26 “ ^^66^ 2F(A^^A22 + Agg 

“ E-l ) - 9(A22B-]5 ■*- -^11^26^ *'* '^66*' -^1 1^22 4.22B-|>|) 


(3.18) 


(3.19) 


Ar 


2B^g(3E^ “ 5Agg) + 2A22(AggF - 3B^gB26) 


^66 


22 '^"66^ 


An •*“ 


+ B22(Ai-jA22 + Agg ~ ) 

^22^^66^22 “ ®26^ 


(5.20) 

(3.21) 


Thus, the governing Ecq. (5.14) for an antisymmetric 
angle-ply laminate is reduced to 

8^ 

ax^Oy^ ■ ^ 3y 




V.' X 


A- ^ + A, — + A5 'T + A, 


7 jjh: = 


(5.22) 
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The displacement components for this type of plate in terms of 
the displacement function, 0 can he obtained from Eqs. (3.10- 
3.12) and are given by 


u 




il£ 






V = 


w 




ll£ 


ay 

8-a r— ^ + a 




5 


+ a. 




g'x -c 'y 

d'^ei 

® 3x^ 3v^ 


'>•2 6 5 i „4 


+ a, 


3y 

3% 

9 3v+ 


v/here 


a .1 — 


ar 


a. 


Biet^Agg - E,) 

® 26^'^66 “ ^® 1 ^ **' ^® 16'^^22 


®26^^22 


a, 


ar 


a, 


®16‘^1 1 




a. 


7 

ag 


Ai ^Agg 

^ 11^22 *^66 " 


ar 


^22^66 


(3.23) 

(3.24) 

(3.25) 


The force and moment resultants in terms of the 
displacement function, 0 can be obtained from the constitutive 
relation (3.2c) and Eqs. (3.23-3.25). These are given by 
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N, 


= a 


^ - 


ifl 




ay 

I a ^0 

y "I 3 3 ,, 33 y 3 

i!£_*„ -ii£ 


H„ = a 


3*^0 


“xy= ^ 


U A 4 2 ■’■ ^15 2 A 4 

ux^r^v ex ov* 


f X 


X oy 


a a s ^0 

^16 '=■ ^-17 ^3i72’ ^18 o‘T'o-T 

' X ' n X n y rj X o y ^ 


19 Sy6 


+ a 


111 


jii 


'V ' ^20 + ^21 * '>22 3/3y4 


«« 


S ^0 

■" ^23 ^ 


y “ ^24 s^5 3y " “25 3^3 9^3 " “26 3^ 3^5 


ill 


+ a, 


d^0 . , 2)60 


(3.26) 

(5.27) 

(5.28) 

(3.29) 

(3.30) 

(3.31) 


where 


^•10 ~ 1-^22 "■ ^“^12^1 ■^12-^66 " ^^66 
+ '^26^^^^^66 ■*■ ^^12 “ 


a^ -j 


^^12^66 “ -^12^1^ 

16^^22^66 

ai2 

“ B^g(3Ai2A66 + 

^1 1-^22 -^12^1 ^ 

“■ ^^26^11^66 

a.j j 

= B^gA22(3A^2 “ 

■*■ -^66^ 



+ B 2 g(- 3 A^ 2®1 **■ ^ 12^66 ^ -^ 11^22 " ^-^66 ''' ^ 
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^14 

- 

® 16*‘^'^66 + "• ■^ 11-^22 ^ 

+ 2 B 2 gA^^Agg 

^15 

= 

^^^ 6 ^ 22 ^ 6 e ■’■ ® 26^®1 ^^66 “ 

22 “ 

^16 

= 

■^ 11^^16 “ ^ 66 ^ 11 ^ 


a ^7 

= 

^^ 11 ^ 16^26 “ ^^ 16-^12 - ^11 ^-^ 11-^22 

. A^g » e 2 ) 



- D ^ 2 -^^^Agg 


^18 


® 16 ^^® 16^22 ■“ ^^ 26 ^ 12 ^ “ ^^^^ 22^66 



" ^ 12 ^'^ 11'^22 ^66 “ ®1 ^ 


a-j g 

= 

^ 22 ^® 16®26 “ ^ 12 -^ 66 ^ 


ro 

o 

“ 

“^ 11 ^^ 16^26 “ -^ 66 ^ 12 ^ 


^21 


■^■ n ^^®26 “ -^ 66 ^ 22 ^ “ 16 ^ 26-^12 “ 

^ 12 ^^ 11^22 





^22 

n 

•^ 22 ^^® 16®26 " ® 12 -^ 66 ^ “ ^® 26-^12 " 

® 22 *‘^ 1 1*^22 



+ Agg - l !*^) 


^23 

= 

•^ 22^^26 ~ ^ 22 ^ 66 ^ 


^24 

- 

^ 11 ^® 16^26 " ^^ 66 ^ 66 ^ ® l 6^^-^66 

~ A ^ g ) 

®'25 

=:= 

3 ( B ^ gA 22 + ® 26 ^ 11 ^ 16^26 ^"^66 

- 3 E ^) 



- ^ 66 ‘■ h ^^ 2 ^ * ^66 - E ?) 


^26 

== 

® 16 ® 26'*'22 " ® 26^®1 ■ ^^ 66 ^ ” ®®66 

^ 22^66 
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3.3.2 Antisyimaetrio Cross-Ply Plates 


It has been shown in Chapter 2 that for an ant isynmie trie 
cross-ply plate 


^16 ~ -^26 " ^16 ~ ^26 ~ ®12 “ ®16 “ ®26 " ®66 ~ ° 
3 Xid . ^22 ~ "*"^11 


(3.32) 

(3.33) 


This GO efficients of the governing Eq. (3.14) for 
this type of plate are reduced, to 


^2 

- -^4 - -^6 ~ -^8 " ^ 


(3.34) 


^ ^66^-^11^11 " ®11 ^ 


(3.35) 

^3 

” ^■^11'^‘’66^ ” "^22^11 ®11^'^11'^22 

+ 4e - E?) 

(3.36) 

A5 

= 2E(A.,.,A22 + Agg E^) - 2E.,B^., 




^66^-^11^22 -^22^11^ 


(3.37) 

Aj 

= 2A22AggE - A.,.,B^,^ + D22(A.j.lA22 

+ Ale ■■ 

(3.38) 

A9 

= Agg(A22D22 “ ®i 1 ^ 


(3.39) 


The displacement components in terms of 0 for anti- 
symmetric cross-ply plates can be obtained from Eqs. (3*10-3.12) 
and are given by 






3x^ 3y^ 





3% 

ax3y+ 


( 3 . 4 - 0 ) 
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i, ^ + K aV + b §1^ 

®7«n 2'\ 2^°Q'\ 4. 

' «x ox C'Y ^ OY^ 


where 


(3.41) 

(3.42) 


^1 ~ ® 11 '^ 66 ^ ^2 ~ ® 11 '^ 22 * ^3 " ® 11 ®r ^4 " ~ ^ 11^1 

h^ = - ^6 ® 11 -^ 66 » ^7 "" ^ 11^66 

bg = 1-422 ■*■ 455 - E-jj Eg = ^ 22^66 


The force and momenf resultants in terms of 0 can be 
obtained from Eq, (3,2c) and Eqs. (3.40-3.42), These are given 
by 



(3.43) 

(3.44) 

(3.45) 

(3.46) 

(3.47) 
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, 3 ^^ , d ^0 , 

27 ?,5 Sy ' ^28 3^3 3 y 3 ■ >’29 ^^ayS 

where 




= 

■^ 12 ^ 66^11 

^11 


B ^ l ( Ai^Ei - - ^ 22 ^ 66 ^ 

^12 

= 

“ -^ 1 1 ^ 66^1 1 

^13 

= 

^ 12 ^ 16 ®!! 

^14 

= 

®11 ^-^ 11^66 " -^ 22 ^ 66 ^ 


~ 

Bii ( Ai 2 Ei + ^6 - 


= 

“ ^ 12 -^ 66^11 

b ^7 


® 1 lV 6^^22 " 

^18 

“ 

^ 12 ^^ 66^11 

^19 

= 


CT * 

ro 

o 


® 11^22 ~ ® 11 ^^ 11^2 ''■ ^66 “ ® i ) “ ■'^11 ■'^ 66^1 2 

^21 

= 

Bii®i " -^ 22 ^ 66^11 ■’ ^ 12 ^^ 11^22 “^66 ~ ^ i ) 

^22 

= 

“ -^ 22 ^ 66^12 

^23 

= 

“ -^ 11 ^ 66^12 

^24 


^ 11^1 ~ ^^ 2^^^^^22 ^66 “ “ ^ 22 ^^^^ 6 S 

^25 

- 

“ Bi 2 '^ 22‘^66 “ ^ 22^^'\'\^22 *** ^66 ~‘ ^ 

^26 

- 

2 

•^ 66^^11 ■" -^ 22 ^ 22 ^ 

^27 

n : 

“ 2^11 


( 3 . 48 ) 



6 ^ 


^28 " " ^^66 ^-^ 11 -^22 
^29 "" " ^^ 66 ^ 22^^66 

In the chapters to follow, eqpiation (3.22) is applied 
to buckling and static problems. 

3.4 METHOD OF SOLUTION 

As stated in Chapter 1, all the problems considered 
hereafter will haye two opposite edges simply supported. 
Consequently , it will be proper to attempt a solution which 
will satisfy these boundary'- conditions identically. The 
displacement function can, therefore, be written as 

0(x, y) = 0^iy) sin(aj^x) (3.49) 

where 

= mVa 

When this expression is substituted in Eq. (3.22), the resulting 
equation for will be an ordinary differential equation and 

will haye a solution of the type 

A OC y 

0^iy) = e “ (3.50)- 

The solution of Eq, (3.22) can, therefore, be written in the form 

0^{x, y) = e ^ sin(aj^x) (3.51) 

Substituting relation (3.51) into the homogeneous part of Eq, 
(3,22), we get 
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A® + + C 2 '/ + = 0 ( 3 . 52 ) 

where 

= Ar^r/Ag j ^2 “ 

Cj = Aj/Ag ; C^ = A^/Ag 

Setting 

= Y 

Eq. ( 3 . 52 ) becomes 

+ Cg + C 3 Y + = 0 (3.53) 

The above equation is a quadratic in Y and can be solved in 
closed form by the well-known Perrari's method. 

Assuming that the possibility of multiple roots is 
remote, the following combination of roots of Eq. (3.53) are 
considered. 

(a) All roots real: 

Let Y = Y^, Y 2 » Y^ and Y^ be the four real roots of 
Eq. (3.53). Then 



= ± Vy^i 

± ^^^2? 

4* 

Vy, 

and 

± % 

or 

A = + A 

+ A 2» 


A ^ 

and 

+ ^ 4 


If all the roots of the characteristic Eq. (3.53) are distinct, 
then eight linearly independent solutions are possible. Thus 
we can write 
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, h V 

^jQ(y) = + do© 


- 


^2 ^ ^ 

2 ^^ d 26 + djS 


jay - y -X.,ay ~>.av 

+ djO 1 + dge 2 ^ 2 “ + dge 


^^3 ^ ^^4 *^117 

4 ^ 


or 


^l“my>+ \ sinh( >.,a^y) 

+ oosh( A 2 V) + ^2 V) 

+ oosh( Sj^y) + sinh( AjCt^y) 

+ oosh( X ^ + \ ^ 4 %y> < 5 . 54 ) 

(b) Tv/o real and tv/o complex; 

Let 

^ = ± VY -,5 + VY2? ± V (Y 5 + i Y4); + Vrf 3""- ' i Y 4 ) 

or A= + A .|5 + },2»i(^3‘*'i^/}.)li(^3“d^4) 

Then the function 0^{y) can be written as 

V) + % >1 V) 

+ oosh( ?12 V) * V ®“h( Aj 

cosh(A3%y) cos(?>^a^y) 

+ cosh(A3a^y) sin(>4ajny) 

+ sinhC A^a^y) sin( A^aj^y) (3,55) 
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(c) /ill roots complex; 

The case of all roots in complex conjugate pairs 
follov/s directly from case (b) and the solution can be 

v/ritten as 

+ sinhC a^y) cos( >^2 

+ Ajjj cosh( A^ttjjjy) sin( 

3 

+ sinh( /-jajiiy) sin( >2“m^^ 

+ V oosbO jC<^) cos( \^a^) 

+ \ slnh(>.,(xj) cos( > 4 a„y) 

+ cosh( ■A 3 %y) sin( \ ^\j) 

+ sinh( A 3 aj^y) sin( A^^jj^y) (3.56) 

The particular solution of Eq, (3.22) depends upon the nature 
of the loading function q(x, y). 

In subsequent chapters, this method is applied for 
the solution of static and buckling problems of layered plates. 


- o 
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CHAPTER 4 

STATIC ANALYSIS OP ANTISWiETRiC ANGLE-PLY PLATES 

4 , 1 INTRODUCTION 

As pointed out in earlier sections, lamination 
asymmetry of geometrical and material properties in laminated 
plates can result in coupling between bending and extension of 
the laminate. This phenomenon is demonstrated by bending of a 
laminate that is subjected to only in-plane forces or extension 
of a laminate that is bent by application of moments only. 

Anisotropic composite plates that are simply supported 
on all sides have been studied extensively . ^ Sun^*^ 
employed the double Fourier series representation for dis- 
placement of an anisotropic plate simply supported on all sides. 
Results for the maximum deflection and bending moments show 

that bending-twisting coupling decreases the plate stiffness. 

61 

Kan and Ito employed the Boussinesq, method to 
obtain integral equilibrium equations of the plate. Single 
Fourier series used in load and displacements reduces the 
system cf partial differential equations to a system of 
ordinary differential equations. In conclusion, authors list 
plate characteristics which influence degree of coupling and 
estimates of the magnitude of the coupling. 

Whitney^^ presented an analysis for bending 
deflections and buckling loads of a simply supported, square, 



67 


single-ply composite plates employing both the Fourier analysis 

and Rita method. He observed that the Fourier analysis provides 

a faster convergence than the Ritz method. 

From the available literature, it appears that 

practically no work has been done to analyse the unsymmetric 

angle-ply plates v/ith two opposite edges simply supported and 

the remaining edges free. In this chapter, the static analysis 

of unsymmetric angle -ply plate is investigated v/ith the above 

boundary conditions. A particular case of this problem is 

presented by Keays who analysed an orthotropic plate with two 

opposite edges simply supported and the other two free. The 

loading, plate geometry and material properties employed for 

solving this problem are assumed to be the same as used by 
63 

Keays with a view to compare the results obtained by this 

63 

technique with that of Keays . The formulation itself is 
quite general and can be applied for the analysis of angle-ply 
laminated plates . 

4.2 ANALYSIS 

4.2.1 Governing Equations 

Consider the case of an antisymmetric angle-ply plate 
of length a, width b and thickness t with two opposite edges 
simply supported and the remaining two free (Figure 4.1). The 
applied load is considered as two line distributions of uniform 
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intensity across the plate width applied at a fixed distance 
from either end of the plate as shown in Figure 4,2. 

The governing equilibrium equation for a thin angle- 
ply plate in bending is derived in Chapter 3 and is given by 






9% 


- 6^ 2 
dx ^y 


4* Ar 




rJx 


■4 Ar 


^ + A ^ 
5x2 9 ayS 


3^J_ 


= aU, y) 


( 3 . 22 ) 


\ 7 here , A^, A^, Aj and A^ are defined in Chapter 3 by 
relations (3.17-3.21). The in-plane forces and moments for 
this type of plate in terms of 0 have been derived in Chapter 3 
and are given by relations (3.26-3.31 )• 


4,2.2 Boundary Conditions 

(a) The boimdary conditions at the ends of the plate are 
those of simply supported edges, namely 

(i) Ih-plane boundary conditions; 

u° = L = 0; X = 0, a (4.1) 

(ii) Transverse boundary conditions; 

w = li^ = 05 x=:0, a (4.2) 

(b) The boimdary conditions at the free edges are given by 

(i) In-plane boundary conditions; 

Hy = F.^y = Oj y = ± b/2 (4.3) 
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(ii) Transverse boimdary conditions; 

?) M a M 

y = ± b/2 (4.4) 

4.2.3 Loading 

The particular load case considered here is a pair of 
line loads syimnetrically placed one inch from either end as 
shown in Figure 4.2. The line loads are assumed to be of 
uniform intensity across the v/idth of the plate. The loading 
being symmetric with respect to the middle of the plate can be 
represented b 3 r a single Fourier sine series containing odd terms 
as 


1 


2 )' ■ 

ra= 1 , 3,5 




where 


= I f sin(a_x) dx 

m 6 m 

For the present case 

Q-m = 9^ ^ mit c/a 

m ^o a 


where 


(4.5) 


(4.6a) 


(4.6b/ 


c = Distance between the line of action of the applied 
load and either end of the plate =1 in, 


4.2.4 Method of Solution 

A solution satisfying the boundary conditions at 
X = 0, a can be written in the form 
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^(x, y) = r ■ sin(a x) (4.7) 

m=1,3,5... 

The solution of homogeneous part of Eq, (3.22) has been 
presented in Section 3.4 of Chapter 3. 

The particular solution can be obtained as follows? 
substituting a single term from the series given by Eq. (4.7) 
into Eq, (3.22) gives 


where superscripts on 0 indicate the order of the derivative 
v/ith respect to y. Substitution of 





'Nq 




into Eq, (4.8), leads to 


(4.9) 


Substituting now the value of from Eq, (4.6b) into Eq, (4.9) 
yields 





' %> 


^ Sin 
8 m 


( 4 . 10 ) 


Incorporating the homogeneous solution obtained in Section 3.4, 
the total solution can be written as 
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0(x, y) = 2 [H^ + A cosh( 

m=1 ,3,5... 1 

+ A alnhC ?, ,%y) + A oosh{ ,X gCi^y) 

^ 5 

+ V slnh(>2«^) + ooshO^jO:^) cosO.^oc,^) 

4 5 

+ V sinhO^a^y) cosCX^ajj^y) 

6 ^ -r 

+ sinhO^ajj^y) sin(>^a^y) 

+ coshO-.^^y) sin(>.. ^“j^y)] sin(a^x) (4.11) 

assuming two roots (X ^ , X. 2 ) ^real and two in complex conjugate 
pair ( /\ ^ ± i X ^ ) . 

Since loading and boundary conditions are symmetric 
with respect to x and y axes, the displacement function J2((x, y) 
will also be symmetric , The asymmetric part of the solution 
is zero and the total solution can then be v/ritten as 


0[x, y) = S. [H„ + V oosh(>,aj^y) + A coBhCX.^a^) 

m — » ^ 

+ Aj^ coshO, a^j^y) cosO . Ojj^y) 

5 

+ sinhCX ^ ctjjpr) sinCX^ajj/)] sin(ajjjX) (4.12) 

Because of symmetiry it is sufficient to satisfy the boundary 
conditions at y = + b/2. The four boundary conditions (4.3) 
and (4.4) give four simultaneous equations v/hich can be solved 
to give \i » \ 5 The in-plane forces and moments 

can be evaluated from relations (3.26-3.31). 
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4.3 COIylPUTATIOUS 

^ *7 

The plate dimensions assumed are^^ 
a = 6.0 in. 

b = 0.75 in. 

t = 0.1 in. 

Typical properties for unidirectional carbon-fibre reinforced 
plastic (GFRP) are 

= 19 X 10^ psi 5 Eg = 3 X 10*^ psi 

G -|2 = 2 X 10^ psi| "^^2 = 0.20 

^21 = '‘^12 = 0*0315 

The stiffness matrices and 3)^^^ for the antisymmetri- 

cally laminated angle-ply plate for various fibre orientations 
are computed by the procedure set out in Chapter 2. The load 
intensity, is set to unity. A computer programme is 
written to evaluate the force and moment resultants at various 
points of the plate. A check on convergence showed that 
deflection of the plate, force and moment resultants converged 
to four figure accuracy v/ithin about 12 terms. On this basis, 
thirteen term solution was used throughout this study. 
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4.4 DISCUSSION AJ4D CONCLUSIO]^S 

The nondimensionalized values of static deflection, w 

and moment M along the centre line of the plate are given in 

Tables 4.1 and 4.2 respectively. The nondimensionalized values 

of moment resultants and transverse shear along a free edge of 

the plate are given in Tables 4.3 and 4*4 respectively. 

However, for the purpose of comparison with the available 
63 

literature , only the dimensional values of v/, M , M and Q 

j y 

are plotted in Pigures 4.3, 4.4, 4.5 and 4.6 respectively. 

The centre line deflection, vf is plotted in Figure 4.3 

for three different fibre orientations (0°, 15° and 30 °) with 

number of layers as a parameter. The figure also shov/s the 

63 

forty-eight term ortho tropic solution given by Keays . It 

can be seen that the plate deflections increase with increase 

in fibre orientations. This may be attributed to the general 

decrease of overall stiffness of the plate v/ith fibre 

orientation for a plate with free edges. The orthotropic 

solution merges into the 0° fibre orientation solution which 

is independent of the number of layers. Further, it may be 

noted that 0° fibre oriented angle-ply plate is analogous to 

63 

the ortho tropic plate of Keays . It is also observed that the 
deflection, w tends to be independent of the number of layers 
when the same is increased beyond N = 4. The difference in 
maximum deflection between four and six-layered configuration, 
for example, is of the order of two percent. This corroborates 
the results of Whitney . ^ ^ 
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The variation of moment resultant, M along the centre 

line and that of M along the free edge of the plate are shown 

in Figures 4.4 and 4.5 respectively for three different fibre 

orientations (0°, 15° and 30°). Again the orthotropic moment 
6 3 

resultants shown in the same figures merge into 0° fibre 
orientation solution. However, the moment resultants unlike 
the plate deflection, vf are found to increase with the increase 
in the number of layers for a given fibre orientation and tend 
to become independent of the number of layers beyond H = 4. 

This may be attributed to the fact that for a given fibre 
orientation, the effective stiffness of the plate and hence 
moment resultants increase v/ith the increase in number of 
layers in a given thickness of the plate. 

A study of the Tables 4. 2-4. 4 show that the moment 
resultants and transverse shear Increase with fibre orientations 
upto 0 = 45° v/hereafter the trend is reversed. This may be due 
to the coupling betv/een bending and extension v/hioh increases 
continuously upto © = 45° and falls progressively thereafter. 

Figure 4.6 shows the variation of transverse shear 
along a free edge of the plate. The trend of variation with the 
number of layers and fibre orientations is the same as for the 
moment resultants . The sudden increase in the value of the 
transverse shear at the point of application of the applied 


load is as expected. 
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Pigures 4«3 to 4*6 shov/ that the orthotropic 
63 

resultants v/ith 48 terms match closely vdth those of 0° 
fibre oriented plate (or anologous ortho tropic plate) with 
13 terms in the series. The convergence of the present 
formulation , is , therefore , faster . 


o 
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TABLE 4.1. VAEIATIOl'T OP NON-DIIffiNSlOFAL DISPLACEIffil'IT (v^EgtV 

ALONG THE CENTEE LOTE OP AN ANTISIMEffiTRICALLY 
LAiroiATED ANGLE-PLY PLATE WITH a/b = 8.0 


Listance 
Along the 
Centre Line, 
X (in.) 

N 

G =15° 
(10"2) 

e =30° 

(10“^) 

0=450 

(10"^) 

0=60° 

(10"^) 

6=75° 

(10“^) 


2 

0.9968 

1 ,6158 

1.9532 

2.8197 

4.1899 

0.50 

4 

0.8638 

1.2459 

1 .2791 

2.2490 

4.0779 


6 

0.8432 

1.1994 

1.2013 

2-. 1788 

4.0609 


2 

1.8925 

3.0679 

3.7044 

5.3513 

7.9540 

1.0 

4 

1.6393 

2.3666 

2.4244 

4.2676 

7.9540 


6 

1.6005 

2.2784 

2.2755 

4.1330 

7.7054 


2 

2.6031 

3.9358 

5.0821 

7.3500 

10.9321 

1.50 

4 

2.2556 

3.2569 

3.3175 

5.8576 

10.6382 


6 

2.2025 

3.1365 

3.1128 

5.6727 

10.6003 


2 

3.1109 

5.0432 

6.0615 

8.7751 

13.0655 

2.0 

4 

2.6965 

3.8950 

3.9490 

6.9897 

12.7146 


6 

2.6330 

3.7509 

3.7034 

6.7678 

12.6673 


2 

3.4161 

5.5381 

6.6484 

9.6237 

14.3455 

2.50 

4 

2.9611 

4.2780 

4.3273 

7.6686 

13.9567 


6 

2.8909 

4.1206 

4.0573 

7.4245 

13.8999 


2 

3.5092 

5.7031 

6.8447 

9.9081 

14.7721 

3.0 

4 

3.0492 

4.4060 

4.4525 

7.8952 

14.3739 


6 

2.9772 

4.2439 

4.1747 

7.6440 

14.3170 
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TABLE 4.2. VARIATION OF NON-DIMENSIONAL MOMENT RESULTANT 

ALONG THE CENTRE LINE OE AN 
Al'ITISYT-METRICALLY LAMHIATED ELATE WITH a/b = 8.0 


Distance 

Along the 
Centre Line, 

X (in.) 

N 

0=15° 

(10“^) 

0 =30° 

(io“^) 

.’0=45° 

(10"^) 

0=60° 

(10-2) 

0=75° 

(10“^) 


2 

0.0946 

0.3406 

0.7941 

0.7992 

0,2308 

0.50 

4 

0.1428 

0.6699 

1.7435 

1.7498 

0,4183 


6 

0.1491 

0.7180 

1 .8854 

1.8649 

0.3933 


2 

0.2989 

1 .0292 

2 . 6867 

4.0383 

2,8916 

1 .0 

4 

0.4168 

1 .7932 

3.1433 

7.6640 

4.3289 


6 

0.4326 

1 .8986 

5.4919 

8.1714 

4.4724 


2 

0.1007 

0.3664 

0.8355 

0.7898 

0.1954 

1.50 

4 

0.1583 

0.7649 

1.9359 

1.8528 

0.4574 


6 

0.1629 

0.8117 

2.0589 

1.8779 

0.34108 


2 

0.0264 

0.1064 

0.2005 

0.0791 

0.0546 

2.0 

4 

0.0523 

0.2911 

0.6451 

0.4213 

0.0980 


6 

0.0553 

0.2977 

0.6455 

0.4224 

0.07158 


2 

0.0056 

0.0271 

0.0320 

0.0445 

0,0756 

2.50 

4 

0.0195 

0.1245 

0,2481 

0.1568 

0.0840 


6 

0.0199 

0.1250 

0.2485 

0.1570 

0.1125 


2 

0.0012 

0.0097 

0.0017 

0.0020 

0.0784 

3.0 

4 

0.0128 

0.0098 

0.1605 

0.1181 

0.0867 


6 

0.0129 

0.0099 

0.1610 

0.1185 

0.1172 



78 


TABLE 4.3. VARIATION OP NON-DIMSIONAL TV/ISTING MOMT 

RESULTANT (M^yE2VA^ ^q^b) ALONG THE FREE EDGE OF AN 
ANTISYMffiTRICALLY LAJHNATED ELATE WITH a/b = 8,0 


Distance 
Along the 
Edge, 

X (in.) 

N 

8=15° 

(-10“^) 

8=30° 

(“10'"^) 

8=45° 

(-lO”^) 

6=* 60° 

(-10-2) 

0=75° 

(-10"2 


2 

0.6631 

2 , 3840 

5.7063 

6.0714 

• 3,8355 

0.50 

4 

0.9669 

4.7308 

11.9332 

13.8008 

5.7103 


6 

1.0220 

4.9612 

13.0458 

15.0476 

6.0569 


2 

0.3799 

1.3729 

3,2410 

3.7910 

2.0039 

1 .0 

4 

0.5647 

2,7207 

6.9363 

7.5636 

2.9877 


6 

0.5988 

2.9675 

7.6203 

8.2623 

3.1701 


2 

0.0932 

0.3459 

0.7421 

0.5962 

0.1693 

1.50 

4 

0.1553 

0.8228 

1.8578 

1.2944 

0.2592 


6 

0.1676 

0.9206 

2.0895 

1.4355 

0.2764 


2 

0.0265 

0.1058 

0.2035 

0.1206 

0.0318 

2.0 

4 

0.0490 

0.2914 

0.5854 

0.2835 

0.0498 


6 

0.0537 

0,3326 

0.6717 

0.3180 

0.0532 


2 

0.0073 

0.0307 

0.0545 

0.0288 

0,0112 

2.50 

4 

0.0144 

0.0944 

0.1732 

0.0687 

0.0174 


6 

0.0160 

0.1095 

0.2018 

0.0772 

0.0185 


2 

0.0 

0.0 

0.0 

0.0 

0,0 

3.0 

4 

0.0 

0.0 

0.0 

0,0 

0.0 


6 

0.0 „ 

0.0 

0.0 

0,0 

0.0 
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TABLE 4.4. VAEIATIOlf OE NON-DIfffiNSIONAL TRANSVERSE SHEAR 

(QyEgtVA^^qQb) ALON& THE FREE EDGE OF AH ANTISYIM- 
TRICALLY LAMINATED AITGLE-PLY PLATE V/ITH a/b = 8.0 


Distance 
Along the 
Edge, 

X (in. ) 

N 

6=15° 

(“10"^) 

©=30° 

(-10“^) 

6=45" 

(-10-2) 

6=60° 

(-10"^) 

6=75° 

O 

(-10~-) 


2 

0.0215 

0.0749 

0.1779 

0.1871 

0.0657 

0.50 

4 

0.0321 

0.1505 

0.3873 

0.3796 

0.0979 


6 

0.0340 

0.1640 

0.4262 

0.4164 

0.1039 


2 

0.1116 

0.4038 

0.9963 

1.3688 

0.8944 

1.0 

4 

0.1552 

0.7113 

1 . 9543 

2.6138 

1,3437 


6 

0.1627 

0.7618 

ro 

• 

O 

2.9082 

1 .4275 


2 

0.0243 

0.0862 

0.2018 

0.2096 

0.0790 

1.50 

4 

0.0369 

0.1770 

0.4463 

0.4273 

0.1183 


6 

0.0393 

0.1833 

0.4924 

0.4690 

0.1256 


2 

0.0087 

0.0325 

0.0711 

0.0642 

0.0286 

2.0 

4 

0.0141 

0.0736 

0. 1710 

0.1577 

0.0443 


6 

0.0152 

0.0818 

0.1911 

0.1519 

0.0473 


2 

0.0047 

0.0180 

0.0396 

0.0434 

0.0285 

2.50 

4 

0.0075 

0.0406 

0.0924 

0.0894 

0.0442 

' 

6 

0.0080 

0.0452 

0.1031 

0,0979 

0.0471 


2 

0.0039 

0.0151 

0.0338 

0.0414 

0.0295 

3.0 

4 

0.0060 

0.0328 

0.0757 

0.0836 

0.0457 


6 

0.0064 

0.0365 

0.0840 

0.0912 

0.0487 








FIG 4 3 VARIATION OF CENTRE-LINE DEFLECTION 







FIG 4 5 VARIATION OF TWISTING MOMENT 
ALONG A FREE EDGE OF THE PLATE 
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CHAPTER 5 

BUCKLING OE ANTISYIvE/IETRIC CROSS- AID AI'IGLB-PLY PLATES 

5.1 INTRODUCTION 

It is well-known that the critical regions of 
aerospace structures utilizing thin plate construction are 
prone to buckling. The application of new design concepts and 
new materials such as fibre reinforced laminated composites 
in aerospace structures have precipitated the need for corres- 
ponding improvements in analysis techniques. A particular 
characteristics of these laminated plates in contrast to 
homogeneous plates is the possible coupling betv/een in-plane 
extension and out-of -plane bending. Such couplings can 
significantly affect the load response characteristics of 
thes6 plates. 

Considerable amount of information is available on 
the buckling of isotropic^^’"'^'* and orthotropic plates 
with various boundary conditions and subjected to various 
types of loadings* However, buckling of laminated composite 
plates has been receiving increased attention in the recent 
past. 

17 

Ashton and Waddoups presented a Rayleigh-Ritz 
solution for buckling of arbitrarily laminated rectangular 
plates subjected to biaxial compression and in -plane shear 
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buckling response of antisymmetrically laminated cross- and 
angle-ply plates for many boundary conditions for which the 
solutions do not seem to be available. 

5.2 BUCKLING MALYSIS 


5.2,1 Derivation of Governing Equation 

For an antisymmetric cross-ply or angle-ply plate 
subjected to an in -plane load N per unit width (Figure 5.1)) 
the governing equation (3.22) can be written as 


3® 


8x 


where 


+ A 


d^0 


= q(x, y) 


d^0 


+ Ar 


aSj. 


+ A, 


9^ 


q(x, y) = “ 

^ c'-x 


( 3 . 22 ) 

(3.1) 


Substituting the expression for w from Eq. (3.25) for anti- 
symmetric angle-ply plate, Eq, (3.22) can be rewritten as 


Ai ^ + A, 
^ c)x^ 


38_£ 


,8 


3 “^ 1 . .2 
oy 


-h A| 






^ ^ <3x^ 3y^ 


+ A, 


9 dyS 


0X^ ^ 3xV3y2 


+ ar 




3x^ Sy"^ 


) 


(5.2) 


A similar governing equation for an antisymmetric cross-ply 
plate can be obtained by substituting the expression for w 
from Eq. (3.42) and is given by 
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A + A -42^ 


+ K 


'd^0 


+ Ar 




+ kr 


hx^ c>y‘- 5 9x^ ^y^ dx^ dy^ 5 ^ 


8 


y 


+ w. 


1::^ + ‘>8 


' ^'x ° dy 


+ br 


3^0 


-) = 0 


(5.3) 


Sx'^ 3y^ 

The coefficients A-| , A^, A^, A^ and A^ in Eq» (5.3) ^or an 
antisyioDietric cross-ply plate are given by relations (3.35“ 
3.39). 


5.2,2 Method of Solution 

The general method of solution is explained in 
Section (3.4) of Chapter 3. The coefficients of the eighth 
order polynomial in X given by Eq. (3.52) take the following 
values for an antisymmetric angle-ply plate; 

C -| = — Arj / Ag 

(^2 ~ '^ 5'^ ^9 ”” 4 g 

C ^ — — A^/ Ag + p.ag// Ag ( 5 . 4 ) 

= A^/Ag — pa^/Ag 
p 

The corresponding coefficients through for an anti- 

syrmietric cross-ply plate are obtained by replacing a^^ 

and ag in the above expressions by by? bg and bg respectively. 
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5.2.3 Boundary Conditions 


The follov/ing sets of boundary conditions are 
considered on the edges y = + b/2 whilst the edges x = 0, a 
are simply supported in each case: 

(i) Both edges free 

(ii) Both edges clamped 

(iii) One edge clamped and the other edge free. 

The boundary conditions for all these cases have been derived 
in tenns of the function 0(x, y) in Chapter 3j» which are 
repeated here for convenience of reference. The coefficients 
of Eq. (3.53) for the buckling problem contain tv^o unknowns 
namely IT and m. The term 'ra' can be identified v/ith the 
number of half waves along x-direction at buckling and take 
integer values depending on the mode of buckling. As indicated 
in Chapter 3j the closed form solution of Eq, (3.53) can be 
readily obtained. Substitution of solution 0j[j(y) given by 
expressions (3.54)> (3.55) or (3.56) in first tv/o of the 
boundary oondltions leads to a set of four elgebraio equations 
in view of the symmetry of the boundary conditions, v/^hile for 
the third set of boundary conditions, a set of eight 
algebraic equations is obtained. The free edge boundary 
conditions, namely 


and 


N. 


N, 


= 0 

= 0 


y 

xy 

SM , 

+ 2 — ^ 

ay 3 X r 


) 

<> 


y 


= b/2 


(5.5) 


° ,/ 
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lead to the following four algebraic equations for an anti~ 
symmetric angle -ply laminate : 


^ ^^ 8 ^Q ~ ^ ( 5 . 6 ) 

^ 2 lVj ^ 22^2 ^ 23 V^ + ° 24 \ ^ 25 ^^^ ^ 26 ^ 

^ 28 Vq ^ ^ ( 5 . 1 ) 

I O 

*^31^ ^ ^52-^m2 ^ 35 \ ^ ^ 34 \ ^ ^35^^ 

+ = 0 (5.8) 

‘^42^112 ^43 ^44^m^ ^45^115 ^46\g 

■*■ ^Al^rj ■*■ ^48Vg == ° 


where 

2 

= P> ^ sinh (S^) [1 - a^23 
0^2 = ^1 cosh (S^) [1 - a^ 2 ^ ] 

0^2 “ ^2 ”" ^123 ^2^ 

C -j ^ “ ^2 ^^2^ "** ^123 ^2^ 

^15 ” sinh (S-^) cos. (S^) ^ ^ 123 ^ 4 ^ 

-• cosh (S^) sin (S^) [ -i- 



cosh (S^) cos (S^) *- 

- sinh (S^) sii (S^) [ 
cosh (S^) cos (S^) [ + 

+ sinh (S^) sin (S^) [ 
sinh (S^) cos (S^) [ + 

+ cosh (S^) sin (S^)[ 

cosh (S^) [1 - > 

sinh (S^) [1 - A 

A 2 cosh (Sg) [l ~ A: 

>2 sinh (S2) Cl “ A! 

cosh (S^) cos (S^) [J^ - 
+ sinh (Sj) sin (S^) 1^2 
sinh (S^) cos (S^) [J-j - 
+ cosh (Sj) sin (S^) [Jg 
cosh (S^) sin (S^) [J^ - 
+ sinh (S^) cos (S^) [- i 


^ 123*^4 ^ 

+ a't23'^3^ 
^123^^3 ^ 

- 23 *^ 4 ^ 

^123^5] 

” ^123^4^ 

h 

;] 

;3 

54*^5^ 

- ai54J6^ 
^154^5^ 

“ ®'154'^6] 

^154^5^ 

2 ^ 154 *^ 6 ^ 
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Ogg = sinh (S3) sto (S^) [J, - 

+ cosh (S^) cos (S^) ["“ Jg ^154^^6^ 

= cosh (S^) [“ 1 + >.-| “• aggo + ^230 ^ 

C^2 = (S-j) [“ 1 + £>'210^1 "" ^"220 ^1 ^230^1 ^ 

= cosh (Sg) [“ 1 + ^2^0^ 2 " ^220^2 ■*■ ^230 ^2^ 

^ 54 sinh (^2^ ^ ^2 10^2 **" ^220 ^ 2 ^2 30 2 ^ 

^55 ” cosh (Sj) cos (S^) [•^ 1 + ^210^1 *" ^220^5 ^250^9^ 

+ sinh (S^) sin (S^) [B^2^0^2 “ ^220^6 ^230^10^ 

C^g — sinh (S^) cos (&^) [*" 1 + ^210^1 ^220^5 ^230*^9^ 

+ cosh (S^) sin (S^) - ^220^6 ^250^10^ 

O^y — cosh (S^) sin (S^) ["^ 1 + ^220^5 ^250^9 ^ 

+ smh (S^) cos (S^) [- a^^QJg + 0^220^6 ^230^10^ 

^38 ^ sinh (S^) sin (S^) 1 + ^210*^1 ^220^5 ^ ^230*^9^ 

+ cosh (S^) cos (S^) C **" ®'210^2 ^220*^6 ^ ^230^10^ 

2 4 (3 

^41 = ^1 ^ ®’240 ^ 1 ’” ^250 ^ 1 ®'260 1 ^ 

246 
^42 “ ^ ’*’ ^240^1 “ ®250 ®'260 ^1^ 
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^43 ^2 ®-240 ^2 ^ ^250 ^2 ''' ^260-^2^ 

*^44 ~ ^2 ^^2^ [- 1 + Ug^Q ^2 ” ^250 '^2 '*' ®’260 ^2^ 

*^45 cosh (S^) sin (S^) *" ^ 250^7 **’ ^260’^11^ 

+ sinh (Sj) cos (S^) [- + ag^QJ^ “ ag^^Jg 

■’' ^260*^12 ^ 

^46 " ^^4^ ^ ^4 ^'240*^3 *" ®'250'^7 ®26o'^11^ 

+ cosh. (S^) cos (S^) ^240^4 ^250*^8 

^260^12^ 

^47 ~ cosh (S^) cos (S^) [ “• ^^ " ag^Qj^ + ®'250*^7 "” ^260‘^11^ 

+ sinh (S^) sin (S^) [“>^ + a240'^4 “ ^250'^8 

^260'^12^ 

^48 ” sinh (S^) cos (S^) ^“^4. *" ^240^3 ^250'^7 ~ ^260'^11^ 

+ cosh (Sj) sin (S^) ^ + ^240^4 "* ^250*^8 

+ ^260*^12^ 

C! _ . q _ S2L 

- 1 2^ » *^2 “ 2 2'R 

q - a5 . S = 

^3 ~ 3 2R ’ 4 

R = a/b = Plate aspect ratio, 


miT 
4 2R 
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^123 

^ !il . 
^12 


^154 

^15 

^14 



®21 . 


^220 

^22 


^210 

^20 


^20 


^230 

= 1?5 „ 


^•240 

(ag^ 


^20 


^^•20 

+ 2a2^) 


( ®'22 

2^26^ . 

^260 

^23 

^250 

Xb.20 + 

^^24^ 

^^20 

+ 2a24} 


Jl = 

Jg = - 2 >. j ^ ^ 

^4 ■“ ^3^ ^3 " ^ '^4^ 

J5 = ^3 +•*! - 6>^>4 

jg = 4:\3A^( - a2) 

J7 = 5 >3 >4(2^3 a| - >|) - 

Jb = 5 A3 >404 - 2 >.3 >4) +^3 

Jg = 15 Xj >, 4 { >.| - > 3 ) + A® - > 4 

J10 = 2 ?-3 V'°^ 3^:4 " ^^3 ■ 

J 11 = 7)3>4(5A5>2 . 3,v^>4 _,.5) ^ ^7 

J 12 = 7 A3X4(5>.3>4 - 3>^ A 4 ->^) + 



96 


The fixed edge boundaiy conditions, namely 

u = v = w = ~|=0 aty = -*b/2 

lead to the following set of four algebraic equations 
ajitisyraraetric angle-ply laminate: 



^52^n2 

^53^112 ■ 

*■ ^54^^ 

" “55\ 

+ °56\ 



^58^ng 

= 0 



°6lV^ 

+ ^62^112 



y 

^66\g 



+ ^68\lg 

:= 0 




‘^72^112 

■*■ ^73^^ ‘ 

'• ^74^4 




■'■ ‘^77‘Vy 

^78\q 

= 0 



°8l\ 

^82\i2 

^83^5 ‘ 

^84^^ 

■*■ ^85V^ 

+ ‘^es^g 



CO 

00 

o 

J- 

= 0 




where 


~ A ■! sinh(S-j ) [“1 **'^21^1 

2 4 

0 ^ ^ cosh (S-| ) f 1 3-2 "j ^ ^51 1 ^ 

2 ' • A 

C = X 2 sinh (>^ 2 ) f ^21 ^ 2 "" ^31 2 ^ 

2 A 

= ^2 oosh (^2^ *” ^2 1^2 ^3 1^2^ 


( 5 . 10 ) 
for an 


(5.11) 


( 5 . 12 ) 


(5.13) 


(5.14) 



0^5 == cosh (S^) sin (S^) [ ] 

+ sinh (S3) cos (S^) [-A3 + - aj-jJg ] 

^56 ~ sinh (S3) sin (S^) [“ A^ ”■ a2-jJ3 + a3-jJY 3 

-1- cosh (S3) cos (S^)[A3 - + a3^jQ] 

^57 ” cosh (S3) cos (S^)CA^ +a2-jJ3 ~a3^J^3 

+ sinh (S3) sin (S^)[ / 3 " a2'|J^ 

G^g = sinh (S3) cos (S^) - a2^J3 + ^31^7] 

+ cosh (S3) sin (S^) ^*“^3 + a2'|<3^^ “ ^ 31 *^ 8 ^ 
^61 " cosh (S-|) [1 “ a^^ A 1 + ag^A^] 

Cg2 = sinh (S^ ) [-1 + **^64^1^ 

^63 ^ cosh. (S2 ) [1 ^54 ^ 2 ^64 ^ 2 ^ 

<^64 - {S2) [~ 1 + ^ 34 ^ 2 ^ ^ 64 ^- 2^ 

0g3 = cosh (S3) cos (S^) [1 “ a^^J^ ■‘■^64*^5^ 

+ sinh (S3) sin (S^) [- ®' 64 '^ 6 ^ 

Cgg = sinh (S3) cos (S^) [-1 + - ag^J^] 

+ cosh (S3) sin (S^) [ag^Jg ” ^64'^6^ 



Og^ = ooah CSj) sin (S^) [- 1 + - ag^Jj] 

+ sinh (S^) cos (S^) [ ^ ^• 54*^2 ^64*^6 ^ 

Ces = sinh (S^) sin (S^) [*- 

+ cosh (S3) COB (S4) - ag^Jg] 

^71 ~ cosh (S-j ) [ 1 " ®'87 ^ 1 ^ ^97 ^ 

Gy2 ~ sinh (S^) [*■ 1 + ^q'j A -j ^97 A -j ] 

Crj,3 = cosh (Sg) [1 ~ agyAg + ag^A^^ 

— sxnh (Sg) f*” ^ ^87^ 2 *” ^97 A 2^ 

^T5 ^ cosh. (Sj). cos ( ) E ^8T^1 ^ 

+ sinh (S^) sin (S^) [- ^ ^97^6^ 

0 ry ^ ^ ^ inh ( 3 ^ ) cos ( ) r ^ 3 q *1 ^97^5^ 

"h cosh s^n ( ) ^^87^2 ^97^6 ^ 

CyY ^ cosh (S^) sin (S^) [*“ 1 + ^87*^1 ^97^5^ 

4* sinh (Sj) cos (S^) E*^ ^87^2 ^ ^97*^6^ 
O^Q ^ sinh (Sj) sin (S^) ^^ r ^97^5] 

+ cosh (S^) cos (S^) [8 .qyJ 2 "^ ^97'^5 ^ 
ag.| = A -1 sinh (S|) C "■ 1 + ag^ A ^ ag^ A ^ ] 
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^82 


^ cosh (S 

r 2 

) [1 "" ^87 

^97 1 ^ 

as 3 


A 2 sinh ( S 

2 ) C *** *1 ^87 ^ 2 

ag?^ ^ ] 

^84 


A 2 cosh (S 

2 ) [ "I "" ag? ^ 2 

ag? > 2 ] 

^83 

r: 

cosh (S^) 

sin (S^)[v\^ + S^QrjJ*^ ““ J ^ 



+ sinh (S 3 

) cos (S^) [-.^^3 

^-87*^4 ■” ®'g7‘^8^ 

^86 

= 

s mh (S 3 ) 

sin (S^) [- “ 

^87^3 ^ ^gY^7 ^ 



+ cosh (S 3 

) cos (S^)[A 3 - 

^87'^4 ^97*^8 ^ 

O 

00 

r: 

cosh (S 3 ) 

COS (S^) f A ^ ^ 87*^3 ^ tj ] 



+ sinh (S 3 

) sin (S 4 ) [ A 3 “ 

®'87'^4 ^ ®'97'^8 ^ 

00 

00 

o 

- 

sinh (S 3 ) 

00 s (S 4 ) E~ A 4 “ 

^87^3 ^ ^97*^7 ^ 



+ cosh (S 3 

) sin (S 4 ) [-A 3 

+ ag?J^ - ®' 97'^8 ^ 

^21 



a3i = 

a^ 

^54 

- 

a^/ a^ ^ 

^■64 

ag/a^ 

as? 

- 

a? ? 

^97 " 

®'g/®'7 


It may be noted that the above eight algebraic equations for 
an antieymmetrio angle-ply plate have been derived assuming 
the most general solution given by relation (3.55) . 
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The corresponding eight algebraic equations for an 
antisymmetric cross-ply plate can be similarly derived by using 
relations given by (5*44-3.45), (5.47-3.48) and (5.40-3.42). 


The determinant of the coefficients of 




given by 


M r 
^11 

ij 

"12 

“13 

"14 

“15 

“16 

0 

— V 

°18 


"22 

“23 

"24 

“25 

“26 

“27 

°28 

;| ^31 

"32 

C33 

"34 

“35 

“36 

C jY 

°38 

i “41 

"42 

“43 

"44 

“45 

“46 

"47 

“48 

*^51 

"52 

“53 

"54 

“55 

“56 

"57 

“58 

"61 

"62 

“63 

"64 

“65 

“66 

"67 

°68 

"71 

"72 

“73 

"74 

“t 5 

“76 

"77 

“78 

! "si 

"82 

“83 

"84 

“85 

“86 

"87 

“88 


set to zero leads to a transcendental, equation in N„ and m. 

The elements of this determinant contain W and m as unknowns. 
For a given integer value of m, the value of which makes 
the above determinant zero is obtained by a method of hit and 
trial. The minimum value of which satisfies this condition 
is the critical buckling load. i. computer programe incorpo- 
rating the various possible solutions given by Eqs. (3*54), 
(3.55) and ( 3 . 56 ) is then written to obtain critical values of 
for antisymmetric oross- and angle -ply plates with various 
aspect ratios, fibre orientations and different number of 
layers. 
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5.3 MULVIERIGAL COMPUTATIONS 

Numerical computations are carried, out for the 
following composite materials: 

(a) Carbon fibre reinforced plastics (CFEP) 

E-i/Er = 6.33f ~ 0*6675 ^12 "" 

(b) Graphite fibre reinforced plastics (GFEp) 

E^j/Eg = 405 = O. 5 O 5 4-^2 = 0.25 

(c) Boron fibre reinforced plastics (BFRP) 

E^l/Eg = IO 5 O-jg/Eg = 0.335 -^-12 = 0.30 

The extensional, coupling and bending stiffnesses have been 
calculated by the method as indicated in Chapter 2 , 

5.4 R ESULTS MB DISCUSSION 

5 . 4.1 Plates With Two Opposite Edges Simply Supported and 
the Remaining Edges Eree 

Non-dimensional buckling loads as a function of 
thickness ratio, t^j/tg of the laminae of a f our-lay ere d anti- 
symmetric angle-ply square plate are plotted in Figure 5.2 for 
three different fibre orientations (30°, 45° and 60°). It is 
seen that for all the orientations, the buckling load reaches 
its maximum at the thickness ratio of 0,4142. This is 
attributed to the fact that the coupling stiffnesses (B^g and 
Bgg) tend to wiish at this thickness ratio (Refer Chapter 2). 
It is also observed that the buckling load is more sensitive 
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to thickness ratios at 9 = 30® than at other orientations. 

This may be due to the existence of maximum coupling effect in 
the region of 0 = 30° (see Table 5.1). The increase in 
buckling load is, therefore, maximum as the critical thickness 
ratio of 0.4142 is approached with a corresponding decrease of 
coupling terms and B 2 g. The variations of non-dimensional 
buckling loads with the thickness ratio for a six-layered 
antisymmetric angle-ply scLuare plate for two different values 
of are shown in Figures 5.3 and 5.4. The buckling load is 
seen to reach its maximum at thickness ratio, R 2 of 2,0 and 
4.1623 at R^=l & 2 respectively for all orientations. The 
coupling stiffnesses B.jg and Bgg have been shown to disappear 
at these thickness ratios for a six-layered laminate (Refer 
Chapter 2). 

The variations of non-dimensional buckling loads as 

a function of aspect ratio for an antisymmetric angle-ply 

plate are shown in Figures 5.5, 5.6 and 5.7 respectively for 

three different fibre orientations ( 30 °, 45° and 60° ) with 

number of layers as a parameter*. It can easily be seen that 

as the number of layers is increased, the buckling load tends 

to be independent of this parameter (N). This corroborates 

1 

the results of Whitney and Leissa for antisymmetric types of 
lay-ups . 

The effect of coupling as a percentage reduction in 
buckling load between four- and two -layered and between six- 
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and four-layered configurations of a square antisymmetric angle- 
ply plate as a function of fibre orientations is shown in 
Table 5 . 1 « The effect is comparatiTely small at extreme 
values of fibre orientations because the coupling stiffnesses 
tend to vanish in these regions. However, the effects of 
coupling are more pronounced in the regions of fibre orienta- 
tions of than those of 45 '^- 75 o. 

The variations of non-dimensional buckling loads 
with fibre orientation, 9 is shown in Figure 5,8 with number 
of layers as a parameter. The buckling loads tend to be 
independent of the number of layers faster in the region of 
450-750 than they do in the region of 150-40° as discussed 
above . 

Figure 5 .9 shows the variation of non-dimensional 
buckling load as a function of aspect ratio for a regular (M = 1 ) 
antisymmetric oross-ply laminate for different number of 
layers in a given thickness. The dependence of buckling load 
progressively decreases with increase in the number of layers. 
Similar observation was made even in the case of antisymmetric 
angle-ply laminates. The variation of buckling load for a 
square antisymmetric cross-ply plate with cross-ply ratio, M 
is shown in Figure 5 . 10 with number of layers, H as a parameter* 
It is seen that the buckling load increases with M, the 
increase being maximum between cross-ply ratios of 1 and 3 
whereafter the buckling load becomes asymptotic with cross-ply 
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ratio. The steep rise in the buckling load in the region of 
m := 1 to 3 could be attributed to the increase in the overall 
thickness of odd’^umbered layers, the fibres of which are 
aligned in the axial direction. Since the overall thickness 
of the plate is constant, the buckling load tends to stabilize 
with increase in odd~nuinbered layers and hence v^ith increase 
in cross-ply ratio. For a designer this may indicate the 
optimal cross-ply ratio beyond which the increase in the ratio 
becomes unnecessary to achieve a given buckling load. 

5*4.2 Plates With Two Opposite Edges Simply Supported and 
the Remaining Edges Clamped 

Figures 5.11 through 5*13 depict the variations of 
buckling load as a function of fibre orientation, 0 for GFRP, 
BPEP and CPRP square plates respectively. It is seen that 
the buckling load reaches a maximum in the region of 9 = 45° 
approximately in each case. It may also be noted here that 
the buckling mode in the x direction changes from 1 to 2 at 
6 > 40° . This feature is distinct from the other two sets of 
boundary conditions where the buckling mode always retains 
the single wave shape. The rising trend for 0° << 9 445° wie-y 
be due to the result of increasing effect of the clamped edge 
moments as 9 increases, thereby increasing the effective 
cross stiffness. This is again distinct from the other two 
cases where such an effect is not present. It may be noted 
that at © == 0° , the magnitude of the buckling load is of the 
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aamG order of magnitude as in the other two cases. This shows 
that the edge conditions are not affecting the huokling load 
as much as at higher orientations. Another point of interest 
is the different magnitude of the buckling load at 0 = 90® 
compared to that at 0=0°. This may be due to the existence 
of a higher mode shape i.e., m = 2 as compared to m = 1 at 
© = 0°. This again is distinct from the other two edge 
conditions. The minima of a two-layered solution in the 
region of © = 20° and © = 70° is also worth noting. A 
composite plate with two loaded edges simply supported and 
the remaining edges fixed is, therefore, at its minimum 
buckling strength for a two-layered lay-up in these regions 
of fibre orientations. As expected, the various solutions 
tend to be independent of the number of layers at extreme 
values of fibre orientations owing to the absence of coupling 
stiffnesses. The results for an infinite number of layers 
are the same as the orthotropic solution in which all coupling 
between bending and extension and Bgg) is ignored. 

Percentage reduction in buokling loads between four- and two- 
layered configurations and those between six- and four-layered 
configurations are shown in Table 5«2 for an antisymmetric 
carbon/epoxy squa;re plate as a function of fibre orientation. 
It is seen that maximum reduction in buokling load due to 
coupling effects occurs at © = 45 ° . However, as the number 
of layers increases from 4 to 6, the buckling loads approach 
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the values for orthotropic plates. This shows that the 
coupling between bending and extension rapidly dies out as 
the number of layers is increased. 


The variation of buckling load as a function of 
thickness ratio, t^/tg for a f our~layered antisymmetric square 
plate for three different fibre orientations (0 = 30°, 45° 
and 60°) is shown in Figure 5.I4. At all orientations, the 
buckling load reaches a maximum at the thickness ratio of 
0.4142 (see Chapter 2), The percentage increase in buckling 
load between an equal thickness ratio laminate (t^/t2 =1) 
and an optimum thickness ratio laminate (t^/t2 = 0.4142) is 
of the order of 8 percent for 30° and 60° fibre oriented 
plates and 10 percent for a 45° oriented plate. This increase 
in buckling load from the design point of view is quite 
significant and this can be treated as the optimum thickness 
ratio for a given design buckling load. Figure 5.15 shows a 
similar variation of buckling load with thiclmess ratio for a 
six-layered antisymmetric square plate for the same fibre 
orientations. The maximum buckling load occurs at thickness 
ratios, Rg = 2,0 and 4.1623 at R^ = 1.0 and 2.0 respectively. 
This may again be attributed to the absence of coupling 
effects occurring at these thickness ratios (see Chapter 2). 

The variations of buckling load with aspect ratio, 

R for two different fibre orientations of 6 = 30° and © = 45° 
are shown in Figures 5*16 and 5.1? respectively for 
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graphite/epoxy ant isyrame -brie angle-ply plates with ntimber of 
layers as a parameter, The cusps in these figures represent 
the changes in the buckling mode shape as the plate aspect 
ratio changes. It can be noted that the curves approach the 
lowest level at each respective curve in the neighbourhood of 
R = 1,0 as the aspect ratio increases. This behaviour is 
typical of an isotropic plate as discussed by Timoshenko and 
Gere^, It can also be observed that the coupling between 
bending and extension has a constant influence of aspect 
ratio for a given number of layers. Figures 5«18 and 5.19 
depict respectively the variations of buckling load of a 
carbon/epoxy composite plate with aspect ratio for two 
different fibre orientations, 9 = 30® and 4-5°. It is 
observed from Figures 5.16 through 5.19 that the pattern of 
cusps is a function of number of layers, angle of fibre 
orientation and degree of anisotropy of the material 
The first change in mode shape pattern, for example, for 
graphite /epoxy plate occurs between aspect ratios 1,25 and 
1,50 for all the number of layers (2, 4' and 6') at 0 = 30°. 

For carbon/epoxy plate, this occurs between aspect ratios 
1.0 and 1,25. However, the location of cusps moves to the 
right from two-layered configuration to a six-layered one. 

The first change in mode shape pattern occurs between aspect 
ratios 1 .50 and 1 ,T5 fer a 45° oriented plate both for the 
graphite /epoxy and carbon/epoxy materials for all the number 
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of layers considered. This may be attributed to the symmetry 
of fibre orientation at 9 = 45°. 

5.4.3 Plates With Two loaded Edges Simply Supported. 

One Fixed and the Remaining Free 

Figures 5.20 through 5.22 show the variations of 
buckling loads for square plates of boron/epoxy, glass/epoxy 
and carbon/epoxy respectively with fibre orientation, 9. It 
is seen that the buckling load progressively decreases with 
increase in fibre orientation. This is in contrast with the 
trend noticed both in the clamped and simply supported edge 
conditions where the maximum load occurs in the neighbourhood 
of 9 = 45°. However, this trend is similar to the case of 
two edges free. The reason, fox this distinct difference can 
be traced to the conditions on the in-plane displacements 
imposed at the edges. In the free edge case there is no 
restriction on the in -plane displacements while these have 
been restricted completely in both the simply supported and 
fixed edge cases. The variations of buckling loads for 
four-layered and six— layered plates are distinctly different 
from those of two— layered plates. However, all the solutions 
tend to be independent of the number of layers at extreme 
values of fibre orientations. The difference between two- 
layered and four-layered solutions as a function of 
cmi be clearly noted. The coupling effects, as expected, 
become more predominant as the ratio increases. For 
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each of the materials considered, however, the coupling 
effects are observed to be maximum in the region of 0 = 45 ° . 
Percentage reduction in buckling load between four and two- 
layered and between six and four-layered configurations of a 
square antisymmetric carbon/epoxy plate is shown in Table 5 . 3 , 
as a function of fibre orientation, 0. A cursory study of 
the Table 5.3 reveals that the coupling effects are more in 
the fibre orientation range of 15 °-' 45 ° than in the correspon- 
ding range of 45 ° '' 75 °. 

The variations of buckling loads with aspect ratio 
for a carbon/epoxy antisymmetric angle-ply plate are shown in 
Pigures 5.23, 5*24 2 nd 5.25 for three different fibre orienta- 
tions, 0 = 30 °, 45° and 60° respectively with number of layers, 

N as a parameter. As discussed earlier, the cusps in these 
figures occur because of changes in the buckle mode shape as 
the plate aspect ratio is changed. However, the aspect ratio at 
which changes in mode shapes occur appear to be dependent upon 
the fibre orientation and number of layers. The first change 
in mode shape for a 30 ° oriented plate with N = 4 and 6 occurs 
between aspect ratios of 2.75 and 3,0 whereas for 45° and 60° 
oriented plates, the corresponding figures are 2.25 and 1.75 
respectively for the same number of layers. Thus with 
increase of fibre orientation, the cusps move backward with 
respect to aspect ratio for four- and six-layered plates. 

The trend for two-layered coupled solution is also similar 
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though numerically different. It is, however, observed that 
for' symmetrical 45° fibre orientation, the cusps become 
independent of the number of layers and occur simultaneously 
at the same aspect ratio. This behaviour which is distinctly 
different from isotropic one, may be attributed to the 
different levels of overall stiffnesses at various fibre 
orientations . 

The effect of thichness ratio, t^/t^ of a four- 
layered ant is5nTime trie carbon/epoxy square plate on the 
buckling load is represented in Figure 5.26 for three different 
fibre orientations, B =30°, 45° and 60® . The buckling load 
is found to be maximum at the thickness ratio of 0.4142 for 
all the orientations considered (see Chapter 2), However, 
the buckling load appears to be more sensitive to thickness 
ratio at 0 = 45° followed by © = 30° and 60° respectively. 

The percentage increase in buckling load between a laminate 
of thickness ratio, ~ ^ identical laminate 

with t^/tg = 0.4142 is of the order of 9 percent for a 45° 
oriented square plate whereas the corresponding figures for 
30° and 60° oriented square plates are of the order of 8 and 
7 percent respectively. An examination of the Table 5*3 
shows that the increase in the buckling load for a given 
fibre orientation is proportional to the coupling effects 
which are more predominant at © = 45° than elsewhere. At the 
critical thickness ratio of 0.4142, since the coupling effects 
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are entirely absent, one can conclude that the Increase in 
buckling load is due to the elimination of coupling. Figure 
5.27 depicts the similar variations of buckling loads for a 
six-layered antisymmetric square plate. The general obser- 
vations with these types of laminates are similar to those 
of a four-layered laminate discussed above. 

Figures 5.28' and 5.29 depict respectively the 
variation of buckling load with cross-ply ratio, M for an 
ant isynraietric ally laminated cross-ply square plates of 
graphite/epoxy and carbon/epoxy materials. The essential 
features are the same as those observed in the case of plates 
with free edges. 

It is now desirable to study the effect of various 
boundary conditions on the variation of buckling load with 
different parameters. Figure 5.30 shows the variation of 
buckling load with fibre orientation for various boundary 
conditions. It is seen that the buckling load falls progre- 
ssively from © s= 0° to © = 90° for all number of layers (2, 

4 and 6) for (a) plates with loaded edges simply supported 
and the remaining edges free and for (b) plates with loaded 
edges simply supported, one side clamped and the remaining 
side free. The fundamental mode shape for each of these 
cases occurs at m = 1. The variation of buckling load, with 
fibre orientation for plates with two loaded edges simply 
supported and the remaining edges clamped is, however, quite 
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different. The buckling loads rise continuously upto the 
fibre orientation of 40° and progressively fall thereafter 
for four- and six -layered plates. For a two-layered plate, 
however, the minimum values occur around © = 20 ° and 70° 
approximately. The distinct change in trend in the variation 
of buckling load with fibre orientation for this type of 
boundary condition may be attributed to the Poisson's effect 
which appears to be predominant at these edge conditions. 

The effect of boundary conditions on the variation of buckling 
load with thickness ratio for a four-layered plate at various 
fibre orientations is shown in Figure 5#31» It can be seen 
easily that the change in buckling load with thickness ratio 
is more sensitive for plates with two loaded edges simply 
supported and the remaining edges clamped than for other types 
of boundary conditions. This effect is, however, found to be 
least for plates with unloaded edges free. 


- o - 
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TABLE 5.1. PERGEiJTAGE EEDUCTIOIT IN BUCKLING LOAD BETWEEN 
DIEEERENT LAKERS OR AN ANTISYMMETRIC CERP 
ANGLE-RLY SQUARE PLATE AS A FUNCTION OF FIBRE 
ORIENTATIONS (FREE EDGES) 


— 1 T 

, , 

Fibre Orientation Percentage Reduction in {Percentage Reduction in 
in Degrees {Buckling Load Between {Buckling Load Betvireen 

}N = 4 and N = 2 |N = 6 and N = 4 

* ' 


15 

6.73 

- 

20 

18.93 

3.24 

25 

22.39 

3.90 

30 

23.92 

4.32 

35 

23.97 

4.20 

40 

22,39 

3.93 

45 

19.67 

3.34 

50 

16.18 

/ 

55 

12.27 

2.23 

60 

8.37 

1.63 

65 

5.36 

1.14 

70 

3.63 

0.98 

75 

1.29 

0.40 
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TABLE 5.2. PERCENTAGE REDACTION IN BUCKLING LOAD BETWEEN 

LIEPERMT LAYERS OP AN ANTISYWIETRIC CPRP ANGLE 
PLY SQUARE PLATE AS A PUNCTION OP FIBRE 
ORIENTATIONS (CLAMPED EDGES) 


, , 

t I , . 

Fibre Orientation |Percentage Reduction in jPercentage Reduction in 
in Degrees {Buckling Load Between {Buckling Load Between 

{N = 4 and N = 2 {N = 6 and N = 4 

I I 

■ - I — - 


15 

9.82 


20 

15.73 

2.20 

25 

18.94 

4 . 04 

30 

25.30 

3.74 

35 

27.03 

5.13 

40 

28.57 

4.80 

45 

30.05 

5.27 

50 

29.31 

5.23 

55 

- 

4.98 

60 

24.63 

4.42 

65 

20.98 

5.77 

70 

16.34 

2.60 

75 

10,33 

2.05 
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TABIilS 5 . 3 . PBROENTACtE REDUCTION IN BUCKLING LOAD BETWEEN ■ 
DIFFERENT LAYERS OP AN ANTISYWETRIC CFRP ANGLE- 
PIiY SQUARE PLATE AS A FUNCTION OF FIBRE 
ORIENTATIONS ( CLAMPED-FREE EDGES) 


• — — 1 

1 

Fibre Orientation Percentage Reduction 

in Degrees iBuckling Load Between 

= 4 and N = 2 

f 

__ L 

t ' 

1 

in '(Percentage Reduction in 
jBuckling Load Between 

JN - 6 and N = 4 

1 

, 1 — . 

15 

12.40 


3.22 

20 

17.57 


3.84 

25 

21 .96 


4.41 

30 

24.78 


4.66 

35 

26.66 


4.79 

40 

27.39 


4.82 

45 ■ 

27.26 


4.50 

50 

26.24 


4.15 

55 ■ 

24.55 

■* 

“ 

60 

21.53 


3.23 

65 

17.76 


2.98 

70 

14.05 


2.12 

75 

7.07 
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CHAPTER 6 

BHCKIiING OE AHTISTTMMETRICALLY laminated , RECTANGtJLAR 
PLATES ON A WINKLER-PASTEBIfAK PODNDATION 

6.1 INTRODUOTTCN 

Till now we have considered the buckling behaviour 

of Icanjnated plates under uniaxial compressive edge loads. In 

this chapter wo propose bo oonsider the behaviour of similar 

plates under both uniaxial and lateral loads dependent of 

transverse displaoemenLs* A similar practical problem of 

great interest is a laminated plate on WinJcler*‘Pastemak 

foundation. The behaviour of a soil performing the support 

function as a foundation has been modelled by Winkler and 

Pastemalc as shown in Pigure 6.1. A considerable amount of 

work has already been done regarding the elastic behaviour of 

isotropic plates resting on elastic foundation. iv is* 

however, only recently that an attempt has been made bo 

combine the coupling eCfoots common with laminated plates in 

44"45 

the presence of an elastic foundation# 

Turvey'^^ presented simple, exact solutions for 
unifornay loaded, simply supported and antisymmetrically 
laminated cross- and angle-ply pla'te resting on Winkler 

Pasternak foundation. In a reoenb pape^:,^ he presented the 
analysis for a uniformly loaded, simply supported, anti 
symmetrically laminated, cross- and angle-ply rectangular 
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plates resting on a similar elastic foundation. The deflections 
and stress couples are evaluated as a function of material 
J coupling and foundation stiffnesses. However, it seems that 

j 

' there IS no study on the buckling response of composite plates 
li resting on V/inkler-Pastemak foundation. The purpose of this 
(chapter is, therefore, to fill this gap. The plates are 
assumed to be simply supported on the two opposite loaded edges 
and vaiious boundary conditions on the remaining two edges are 
considered. Tho influence of bhe foundation, material coupling 
a 2 id boundary conditions on bhe buckling response of anti- 
symmetrically laminated cross- and angle-ply plates is studied. 


6.2 BUOigiHCr MAJjYSIS 


6*2 . 1 Ant isvmme brio Anfd.e-ply Plates 


The governing eighth order partial differential 
equation (3.22) as derived in Chapter 3 con be written as 


A 


1 






3??y^ 


+ 







= q (x, y) 


(3.22) 


For an antisymmotrio angle -ply plate resting on Winkler- 
Pastemak foundation and subjected to an in-plane load per 
unit width as shown in figure 6.1, the loading term q(x, y) 
in the above equation can be written as 
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(i(x, y) = 


N 


o^w 


kw + g 

Ox 


-.2 

+ eji\ 

n 2'' 

'-y 


( 6 . 1 ) 


where 

k = WinJcler foimdation stiffness 
g = Pastemalc foundation stiffness. 

Substituting the expression for w from Eq. (3.25) into Eq. 
(3.22), the resulting eighth order polynomial in is 

-h [ - Ay/Ag - GD^^ag / P^Ag] 

+ [Aj^/Ag ~ M’O'g/Ag + G(ag + ag)D^^/p^Ag+ KD-j ^Sg/P ^Ag] 

+ C“ *" G(aY + ag)D-j^/P Ag ■“ KD^^ag/p A 

■I* [ A-^ / Ag “* p s-y/ Ag + GD ^ ^ Q-y/ P Ag + ICD ^ -j &tj/ P Ag ] = 0 

( 6 . 2 ) 

where 

K = kbVr>i -1 

G = gbVD-11 
p = mit/R 

Equation (6.2) can now be solved by the method as discussed in 
Chapter 3. 


6.2.2 Antisymmetric Cross-plv Plates 

An eighth order poljmomial in A for on antisymmetric 
cross^ply plate corresponding to Eq, (6,2) developed for an 
antisymmetric angle -ply plate is obtained by replacing ag 
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and ag in Eq, (6,2) by by, bg and bg given by relations (3.75“ 
3.77) respectively. The coefficients A^, A^, Ay and Ag 
in Eq. (6.2) are replaced by those of antisymmetric croBS“ply 
plates. These are given by relations (3.35) to (3.39). 

The coefficients of Eq. (6.2) for given values of K 
and G contain two unknowns and m. m corresponds to number 
of half waves that the plate buckles. In fact m is an integer 
fond taJces on the values 1, 2 etc., depending on the mode of 
buckling. The estimation of buckling loads for given boundary 
conditions is then carried out by the procedure outlined in 
Chapter 5. 

The following sets of boundary conditions have been 
considered; 

(a) Two loaded edges simply supported and the other two 
edges free. 

(b) Two loaded edges simply supported and the other two 
edges clamped. 

(c) Two loaded edges simply supported, one edge clamped 
and the remaining edge free. 


6 • 3 RESULTS MP DISCUSSION 

The numerical computations are carried out for the 
following data of carbon fibre reinforced plastics; 

E^/B2 = 6.335 G^2/®2 = 
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The variation of buckling load as a function of nondimensional 
Winkler foundation stiffness K for three different fibre 
orientations (30% 45° and 60% of angle-ply square plates 
with two loaded edges simply supported and the remaining sides 
free are shown in Figures 6 , 2 , 6.3 and 6.4 respectively with 
nondimensional Pasternak foundation stiffness, G and number of 
layers, N as parameters. It is seen that a small change in 
the value of the Pasternak foundation stiffness results in a 
greater chahge in the values of buckling loads than does an 
equal change in the values of the Winkler foundation stiffness. 
The percentage change in the buckling loads for a two-layered 
configuration, for example, between G- = 2.0 and G = 4.0 at 
K = 2,0 is of the order of twenty whereas the corresponding 
figure between K t= 2,0 and IC = 4 .O at G = 0 is only of the 
order of three percent for a 30° fibre oriented plate. The 
order of these changes appears to be fairly independent of the 
fibre orientations. Further the buckling load increases 
linearly with K. Figure 6.5 shows the variation of buckling 
load with Winkler foundation stiffness for an antisymmetrically 
laminated cross^ply square plate with two loaded edges simply 
supported and the remaining sides free with G and N as para- 
meters, It is observed that the buckling load for a two- 
layered laminate at G = 5.0 approaches that for a four^layered 
laminate at G = 3,0, This may be attributed to the general 
increase of overall stiffness of the plate with increases in 
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Kj G and W, Howevor, "the off act of coupling rapidly diminishes 
as the nuiubor of layers in the laminate is increased. 

The variations of buckling load with Winkler 
stiffness j K for an antisyimnetric square plate are shown in 
Figures 6.6} 6»7 and 6,8 for throe different fibre orientations} 
© = 30° , 45° and 60° respectively with loaded edges simply 
supported and the remaining sides clamped. Here the buckling 
load for a 30° oriented, two-layered laminate at Pastemalc 
foundation stiffness values of 6,0 and 8,0 are seen to approach 
those for a six-layered laminate of the corresponding foundation 
stiffness values of 0,0 and 2,0 respectively. However, no such 
overlapping of buckling loads is observed for 45° and 60° 
oriented plates. Thus changes in buckling loads for given 
changes in Pastemalc stiffness values or in Winkler stiffness 
appear to be dependent upon fibre orientation for this type of 
boundary conditions. The percentage change in buckling load, 
for example, for a two-layered laminate between G of 2,0 and 
4.0 at K = 2.0 is of the order of eleven for © = 30°, five for 
© = 45° and 2.5 for © 60°. The corresponding changes in 

buckling loads due to changes in Winkler stiffness values from 
K = 2.0 to K = 4.0 at G = 0 are of the order of 0.55 percent 
;Por © = 30°, 0,10 percent for © = 45° and 0,05 percent for 
0=:6O° , Thus the changes in buckling loads for given changes 
in Winkler or Pastemalc stiffness values are dependent upon 
the fibre orientation for a plate with loaded edges simply 
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supported and the remaining edges clamped, though the earlier 
observation of dominance of Pasternak stiffness in comparison 
to Winkler stiffness remains valid even in this case« further, 
for a given value of G, there is very little change in the 
buckling load with K. 

figures 6,9 through 6,11 show respectively the 
variations of buckling loads with V/inkler stiffness, K for an 
antisymmetric square plate with tv/o loaded edges simply 
supported, one side fixed and the remaining side free for 
three fibre orientations, © = 30^^, 45° and 60° and with number 
of layers (N) and Pasternak shear stiffness (G) as parameters. 
An examination of these figures reveals that an addition to 
the general dominance of Pastemalc stiffness over Winkler 
stiffness (in terms of increase in buckling loads), the 
changes in buckling loads appear to be dependent upon fibre 
orientations. The percentage increase in buckling load for a 
two~layered laminate, for example, between G = 2,0 and 4,0 at 
K = 2 is of the order of fifteen for © = 30°, twelve for 
© = 45° and nine for © = 60°, The corresponding figures 
between K = 2,0 and 4,0 at G = 0 for a two^layered laminate 
axe 1.9, 1,5 and 1.2 percent for © =30°, 45° and 60° 
respectively. 

The variations of buckling load with fibre 
orientation, © is shown in Figure 6,12 for a square anti- 
symmetric angle-ply plate with two opposite edges simply 
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support©d EUid bhe remaining edges fixed and with number of 
layers j Winlcler and Pasteinak foundation stiffnesses as para^ 
meters. It is of interest to note that both foundation 
stiffnesses (K and G) exert greater influence on buolcling load 
in the fibre orientation range of 5° to 25° than at other 
fibre orientations. The variation of buckling load for a two- 
layered laminate at K = 10,0 and G = 5»0 follov/s partly a 
different trend than v/hen the Pasternak foundation is absent. 
However, as expected, the solutions at different foundation 
stiffnesses tend to become independent of the number of layers 
at extreme values of fibre orientations. The variation of 
buckling load with fibre orientation for a square ant is 3 rmme trie 
angle-ply plate with two loaded edges simply supported, one 
side clamped and the remaining side free is shown in Figure 
6,13 with number of layers (H), V/inkler stiffness (M) and 
Pasternalc stiffness (G) as parameters. The general observations 
for those variations appear to be common with those discussed 
above . 

A comparison of Figures 6.2 ~ 6,13 for different 
boundary conditions reveal that the increments in buckling load 
levels for given increments of K, G and N are maximum for 
free-free boundary conditions and minimum for clamped^clamped 
boundary conditions. For other boundary conditions, the 
values lie in between. Furthermore, the effect of fibre 
orientation on the buckling load is quite signif ioant for some 
types of boundary conditions. 
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It can be concluded from the above study that both 
material coupling and foundation stiffnesses affect the buckling 
response of the antisymmetric cross- and angle-ply plates 
significantly. The effect is more significant for plates with 
free edges than for those with fixed edges. The changes in 
buckling loads due to changes in foundation stiffnesses also 
appear to be dependent upon fibre orientations for certain 
typos of boundary conditions. It is further observed that a 
change in the Pasternak foundation stiffness exerts a far 
greater change in tho buckling load than does a similar change 
in the Winkler foundation stiffness. 
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FIG. 6-1; GEOMETRY OF THE PLATE RESTING 
ON WINK L E R- PAST ER NA K F'O UN 0 A,T 1 0 N. 




NoAdtmensiQ-nai buckitn^ load 


M'S 

N 

LI 

o 4 0 

X 

z 


/ -I » 


20 

^^ 123456 / 0 ^ 

_ ; , ^ K — : ’ . ' ' 

FIG 6 2 VARIATION OF BUCkL^NG LOAD WITH WiNKtgJ 
STIFFNESS, K Foi? A SaUARE ANTISYMMETRIC ANGLE-ptY PLATE. 




ORdim«nsionaf 



K 

FIG 6 3 VARIATION OF BUCKLING LOAD 


WITH WINKLER 


STIFFNESS, K FOR A SQUARE ANTISYMMETRIC ANGLE-PLY PLATE' 





“O o 


=: 6 0 


H 5 

pG 6 6 \/ADi K ® 9 10 

lOUWDATlOW NONOIMENSIONAL BUCKLING LOAD WITH \a/imi/» cn 






f'lG 6 7 VARIATION OF NONDIMENSIONAL BUCKLING LOAD WITH 
WINKLER FOUNDATION STIFFNESS FOR AN ANTISYMMETRIC 
REY SQUARE PLATE WITH NUMBER OF LAYERS (N)& PASTERNAK 


stiffness (G) AS PARAMETERS 



G .-10 



1 


!tO 

4 -> 

(M 

UJ 


c. 

o 

X 

2 







*“* J'*<C 

-0 

^ 





jL.<>_^^ 





19 

C O 
O 

c 

m c 

6 c 

•ooc 

^ o 
O - 


^"5 0 


0 

4 






2 x » 


40 


3 0 t 




:= 0 


= 2 
= 0 


J . 

4 


®«lIr Yoto™? io°« "irn ' 



Nondimensrona f 
buckling load. 


( \ ! 



2 50j ——-I 1 1 1 1 1 i 1 I 

1 2 34 56789 10 

K 

FIG 610 VARIATION OF NONDIMENSIONA L BUCKLING 
LOAD WITH WINKLER FOUNDATION STIFFNESS FOR A S8UARE 
ANTISYMMETRIC ANGIE-PLV PLATE WITH NUMBER OF 
LAYERS (N) a PASTERNAK SHEAR STIFFNESS (g) AS PARAMETERS 



Nondimens.onal buckhns load 


2 

Nzz6 


C f {( p 

^-10 
10 0 00 " 


ro 

Ni 

^ /V ) 


z 3 5 


with WINKLER BUCKL/NG LOAD 

antisymmetric ANGtF-Piv ^ SQUARE 

(N) & PASTERnL\hFAR 

PI CKNAK SHEAR STIFFNESS (g) AS PARAMETERS 


c 



, Fibre orientation, 6 «■ 

nG 6 12 VARIATION OF NONOIMENSIONAL BUCKLING LOAD WITH ANGLE 
-Of fibre ORIENTATION FOR A SQUARE ANTISYMMETRIC PLATE WITH 
DUMBER OF LAYERS (N), WINKLER STIFFNESS (K) AND PASTERNAK 
IHEAR stiffness (G) AS PARAMETERS 


Nondimenstona! 
buckling load, 


\i 


rO 

■M 

(M 

IxJ 




5 


4 


3 


2 


0 


10 20 30 40 50 60 70 80 90 


Fibre orientation j 0 

FIG 6 t3 VARIATION OF NONOIMENSIONAL BUCKLING LOAD WITH 
ANGLE OF FIBRE ORIENTATION FOR A SQUARE ANTISYMMETRIC ANGLE- 
PLY PLATE WITH NUMBER OF LAYERS CN), WINKLER STIFFNESS (K) 
AND PASTERNAK SHEAR STIFFNESS (G) AS PARAMETERS 


169 


CHAPTER 7 

BtlOIOiING OP MTISyMMETRIC CROSS- MD ANGLE-PLY 003VIPOSITE 
PLATES SUBJECTED TO VARYING UNI-AXIAD LOADS 

7.1 INTRODUCTION 

The next logical set of studies on buckling of 

laiainabed plates will be to extend the analysis to combined 

and varying edge loads. Study of literature shows that 
13 

Ohamis analyzed Ihe buckling of simply supported aniso- 
tropic plates under combined loading* membrane, normal and 
shear loads. He also investigated the effects of the various 
coupling responses on the buckling load. There seems to bo 
no dnf 01 matron on the buckling of anisotropic laminated 

'I 

, plates subjected to vaa:ying in-plane loads# The puipose of 

,1 

!' this chapter is, therefore, to develop an efficient formu- 

^iation to predict the buckling response of laminated plates 

i j 

Vsubjeotod to varying in-plane loads at the edges. It is 
somewhat difficult to attempt a closed form solution direotly 
from bhe differential eq.uation due to the fact that it is 
difficult to satisfy the boundary conditions on the loaded 
edges. Ono may, however, get a closed form solution, in 
principle, by expanding the edge loads in the form of a 
series and matching it with the series solution of Bq,, (3.22). 
This, however, leads to an infinite set of equations and 
consequently a oharaoteristio matrix of infinite size, Eor 



170 


nuirioriCi?-! computatioiis ono isj therefore , forced to malce a 
choice of the number of terms and limit the size of the 
matrix for a viable solution consistent with the needs of 
Gon.vergonce« Consequently it turns out to be like any other 
approximate teclmiquc. In this study, the well-known G-alerkin 
technique v/as, thorcifore, chosen to study the problem, The 
numerical solution is obtained for a particular case of 


voriablo loading, namely/- pure bending by the vvell-lmovAi 


Power method 


74,75 


The formulation leads to an iterative 


convOD?p;once which approximates the buckling loads as closely;' 
as desirod . 


7 . 2 

7.2,1 iurul symmet rlo Mgle-ply Laminates 

Let a simply supported rect^gular plate as shown 
in figure 7,1 be subjeoted to a system of distributed forces, 
acting in the middle plane of the plate. The distribution of 
these applied forces may be expressed by the relation 

14^ = (1 - a y/b) (7.1 ) 

whore Nq is the intensity of compressive force at the edge 
y = 0 and a is a numerical factor. It can be seen that by 
changing a, various particular cases of practical interest 
ai'e obtained, a = 2, for example, corresponds to the case 
of pure bending. 
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The govoi’ning differential equation for an anti- 
eyra^Tietric angle--ply plate is given by (Chapter 5) 

A„ -V A, -iiX, . A. -^X . A, 




•^^1 ,■ 8 '■ ^3 X 6 ’ V 2 ' ^"5 

‘ X ' X '-y 




+ ”8 = ^ 


( 3 . 22 ) 


For the plate loading shovai in Figure 7.1) i^n. (3.22) becomes 




Ar 




’h Ar- 


.il 


r) 

An ,: — ^ 'I 


., 2 

n V 


Q -77 + IT,, “--v = 

9 :•., X A ^^2 




(7.2) 


Substituting the expression for from Bq. (7.1) and ’<v from 
Eq, (3.25) i^^to Bq. (7.2), we have 


>, 8,-/ 13 8 h 


8, 


*■ ■ *, • -V • *» 5, 


‘ ’X 




(1 - ot yA) 0 + agfoO " “ yM 


J 

ih 

dy' 


-i- ac,ITo(1 - ay/b) ^ 

- ^ '-.'x 8y 


(7.3) 


The present problem is solved for the following set of simply 
supported edge boundary conditions: 

w = 0 , = 0 


o 


u =. 0 , 


n = 0 

xy 


X = Oj a 
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w = 0 , 

= 

0 

1 


0 

y 


y = 0, b 

V = 0 , 

11 

0 

I 


By the choice of 


(7.4) 


0 ~ ,> a sin miix/a sin nixy/b (7.5) 

]n=:1 n=1 

the boundary conditions (7*4) are identically satisfied. 

This con bo verified by the expressions of displacements and 
stress resultants in terns of 0, which are given by relations 
(3.23) through (3.31). 

Now applying the Galerkin technique and substituting 
for 0 from Eq. (7.5) into Eq, (7.3), we get 


t f C 5- 2. Bln ES . eta aS] sin Ef . sto as 

6 0 m n 


a b 


•Ho / f [ 1 r nl. S|S . oln E|2 . 


0 0 m n 


sin dxdy -i- 


I a >- I- B^a^y sill 


0 0 


HI n 


a 


sin sin . sin ^ dxdy = 0 


(7.6) 


p, q = 1 , 2, 


where 


4 4 

m'K\ /nn\ 


= A, (fS) + AjC-) (~) . Ajl— ) ) 

2 6 8 
A /ns\ . A 

+ A.7(— ) (-^— ) + A9I13 ; 



Equation (7.6) yields a set of iroai equations to determine 
Making use of the identity 


a 


sin . sin dx = 0 for i 3 

a a 


= a/2 for i = 0 

the first term of Eq, (7.6) can be written as 


Bi 

r 


b 


f / '2- V 

m o t n 


V - 


sin : sin dy ] . 


pTtx 


sin dx 
a 


.r I [| vv ®r 

0 m ^ ^ 


[ / z W ste . sto ES ax] 


2 i- ~ ‘^;;r “mq^mq 

0 m 


b a 




Applying this simplification to Eq. (7.6), the 


result is 


a b 


f Apq-pq - W ”o I / / t 1. I V • 


(7.7) 
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The lac't terra ia Eq, (7.7) can be further simplified as 
a b 


^ o / *■ 'm n ^ sin^] 


. sin 


RSx 


a 


3 in 


aiLz 


dxdy 


{ 


■S’ 

n 


a 

[ I (.2 

0 


m 

. sin 


a 


b 


o n 


aTty 


sin 


aM 


2 i "pn^pn*' b b 


dy 


a 


n 


B a Tv Sin -if 
•^pn pn ^ b 


llSZ sin a^ (iy 


~ '7' B a I 

2 ■“ pn pn qn 


where I 


qn 


; y sto as sin ^ 

o 


dy 


b 

^ ’ 


q = n 


0, q/aj q+nisan eyen number 
2 


^ - ib 


m- 




Ti (a •" q ) 


• q a, q + a 

is an odd number 
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Thus Eq.. (7.7) can be written as 

f - fo f Bp,-p, 


■" '^0 f f ( f;, 


= 0 


p — 1j 2j «.»• IHy •«•*» 

q. * 1> 2^ •••• Uj ««•• 


(7.8) 


This leads to an infinite set of homogeneous equations for 
^pq* assumed that this has a solution. The corres- 

ponding matrix say A can be solved by a process of iteration 
to determine N^. The determination of eigen values of the 
matrix A is made by the well-known Power method. 


7.3 HUMERiaAIi OOMPUTATlONg 

Numerical computations are carried out for two 
different materials, namely carbon fibre reinforced plastics 
(OPRP) and graphite fibre reinforced plastics (G-PRP) with 
the following datas 

CPRP 

E^/Eg = 6.333? " 0.667? 

G-ERP 

E^/Eg == 40,0? ^ 12/^2 = 


'>>^2 = 0*20 


*^12 ~ 
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A gGiiGralizGd coinputGr programine is v/ritljGn with a provision 
for considGration of different values of a. However) in view 
of the multiplicity of parameters that characterize the 
px^esent problem for a given loading condition) only the case 
of pure bending corresponding to a = 2 is considered for 
antisyimnetric angle-ply and cross-ply plates. 

7 . 4 HBSHLfS MD DIS G USSIOl? 

The variation of buckling load with fibre orientation 

is shown in Figure 7.2 for a square antisymmetric graphite/ 

epoxy ohgle-ply plate with number of layers as a parameter. 

These are reduced to uniaxial uniform compressive case for 

1 

comparison with Whitney and Leissa . It can be seen that the 
coupling tends to x’educe the buckling load. The results for 
an infinite number of layers in this figure are the same as 
the orthotropic solution in which all coupling between 
bonding and extension is ignored. For a square two-layered 
plate, tho actual results are about 65% below the six- 
layered results nt 0 - 45°. A similar phenomenon was also 
observed by Jones for antisymmetric cross-ply square plates 
of the same material (GFRP). It is also observed that as 
the number of layers increases beyond two, the buckling loads 
rapidly approach the orthotropic solution. This shows that 
the coupling between bending and extension rapidly dies out 
as the number of layers ificreases. The percentage difference 
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in "the budding loads between a sixdayered plate and a four~ 
layered one, for example, is only of the order of 10 as 
against 61 between four“ and two-layered square plates at 
© = 45°. Figure 7.3 shows similar results .for a square anti- 
symmetric carbon/epoxy plate. 

Figures 7«4 through 7.6 show the variations of 
buckling load with aspect ratio for a carbon/epoxy anti- 
symmetric angle-ply plate for three different fibre 
orientations of 30°, 45° and 60° respectively. The cusps in 
these figures occur because of the changes in the mode shape 

as the plate aspect ratio changes. It can be noted that the 

occurs 

location of the cuspi^at an aspect ratio which is a function 
of the fibre orientation and number of layers. The first 
change in mode shape for a 30° fibre orientation occurs between 
aspect ratios of 1.10 and 1.20 for a two-layered laminate 
whereas the corresponding change in mode shape occurs between 
aspect ratios of 1,20 and 1.30 for both the four-layered and 
six-layered laminates. However, because of symmetry of fibre 
orientation at & = 45°, the first change in mode shape occurs 
between aspect ratios of 0.90 and 1.0 for all the number of 
layers (H = 2, 4, 6). For a 60° fibre oriented plate, the 
above figures are approximately between 0.70 and 0.80 for 
all the number of layers considered. Regarding the effect of 
Goupling between bending and extension, it can be seen that 
this is constant irrespective of aspect ratio for a given 
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number of layors which means that the curves in these 
figures have the same relation to one another, The variation 
of buckling load with aspect ratio for a regular anti- 
symmetric carbon/epoxy cross-ply plate are shov/n in Figures 
7.7 and 7.8 for two different loadings, namely a = 2 (pure- 
bending) and a = 0 (uniform uniaxial compression). The 
comparison of these two figures reveals another interesting 
feature, namely that the location of cusps besides being 
dependent upon fibre orientation and number of layers (as 
noted in the case of antisymmetric angle-ply plates above), 
also depends upon the in-plane loading. The first change in 
mode shape, for example, for a = 2 occurs between aspect 
ratios of 0.90 and 1.0 whereas the corresponding figures for 
a - 0 lie betv/een 1,4 and 1.5 for all number of layers 
considered. The number of changes in mode shape is also more 
for a: = 2 case than for a = 0 case. 

Figures 7.9 through 7.12 depict the variations of 
buckling loads with thickness ratios of the laminae consti- 
tuting a laminate of a given thickness for two different 
in-plane loading conditions (a =0, 2). The variation of 
buckling load with thiokiess ratio, t^/t 2 is shown in Figure 
7.9 for three different fibre orientations 50°, 45° and 60° 
for the case of pure bending. The percentage increase in 
buckling load for a 45° oriented plate between an equal 
thiokness laminate (t^/t^ = 1.0) and a critical one (t^/tg = 
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0.4 I42) is of tho order of 10.5 whereas for 30° and 60° 
plates, these figures are of the order of 8.7. Thus the use 
of critical thickness ratio laminates from the design point of 
view is quite advantageous, considering the substantial 
increase of buckling loads that such laminates accrue to the 
designers. The maximum advantage for a 45° oriented plate is 
obvious from the maximum coupling effects associated with 
this fibre orientation. The variation of buckling load with 
thickness ratio for a six-layered laminate for the case of 
pure bending is shown in Figure 7 . 10 for three different 
fibre orientations and for = 1.0 and 2,0. It is seen that 
the maximum buckling load for R^ = 1 ,0 curve occurs at Rg “ 

2,0 whereas for R^ = 2.0, the maximum load is obtained at 
R^ = 4.4162 (see Chapter 2). Figures 7»11 and 7.12 shovr the 
variations of buckling load with thickness ratio for the case 
of uniform uniaxial loading (a = 0) for four- and six-layered 
plates respectively at two different fibre orientations of 
30 ° and 45°. A comparison of Figure 7.9 for uniform uniaxial 
oomprossion loading v/ith that of Figure 7*10 for pure bending 
shows that the maximum increase in buckling load in both the 
cases occurs at flbro orientation of 45°. However, the 
percentage increase in the former case is marginally more 
than in the latter case. 

The variations of buckling load with cross-ply 
ratio, M for an antisymmetric square cross-ply plate are 



180 


shomi in Figures 7.13 and 7.14 for two different materials 
namely GPRP and GFEP respectively for the case of uniform 
uniaxial loading. It is observed that the increase in 
buckling load may also depend upon the material and there- 
fore the degree of anisotropy. The percentage increase in 
buckling load, for example, between cross-ply ratios of 1.0 
and 10.0 is of the order of seven for CPEP whereas it is 
14 percent for GFRP for a four -layered plate. 

- o - 
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CHAPTER 8 

EESUIjTS MD GORCLUSIONS 

8 . 1 & EWEm oomusioi^Ts 

Two aspects of composite plates namely (i) character- 
isation of composites and (ii) buckling of composite plates 
with various boundary conditions subjected to uniform uniaxial 
and varying in-plane loads with or v/ithout lateral loads have 
been considered in this thesis. 

An attempt has been made to eliminate the effect of 
coupling between in-plane extension and out-of-plane bending 
of composite plates by a choice of suitable thickness ratios 
of the laminae of a given laminate and thereby reduce the 
degradation in the overall stiffness of the laminate. It is 
shown that for a four-layered antisymnetric angle-ply 
laminate there exists a thickness ratio (t^/t^) of 0.4142 at 
which the coupling between in-plane extension and out-of -plane 
bending is completely eliminated. This unique thiclcness ratio 
is further shown to be independent of fibre orientations and 
material properties. Eor other types of laminates, it is 
shown that the degrading coupling effect can in general be 
alleviated or partly reduced. However, for a particular 
case of stacking sequence of ( 71/4 + ©), it is shorn that 
there exists a thickness ratio for a six-layered laminate at 
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which the coupling is entirely eliminated. On the other hand, 
for six-layered antisyraraetrio angle-ply laminates, it is shown 
that there exist infinite combinations of (= t^/t^) and 
Rp (= tp/tj) which lead to special orthotropic conditions 
(complete elimination of coupling stiffnesses). 

The static and buckling analyses of composite plates 
are based upon the solution of governing eighth order partial 
diff ei'ential equation which is obtained from the three 
governing equations of composite plates. Both antisymmetric, 
cross- and angle-ply plates are considered* In order to 
confirm the validity of the proposed analysis, a specimen 
problem on the static response of antisymmetric angle-ply 
plate subjected to transverse line loads is solved and the 
results are compared with those existing in literature. 
Moments and reactions are found to increase with increase in 
fibre orientation upto © - 45® whereafter the trend is 
reversed. This is attributed to the coupling between bending 
and extension which increases continuously from 9 = 0° to 
9 = 45° and decreases progressively from 9 = 45° to 9 = 90°. 

The variations of buckling loads of antisymmetric 
angle— ply plates with fibre orientations for various boundary 
conditions reveal an interesting features. The buckling 
loads decrease continuously from © = 0° to © = 90° for plates 
with two loaded edges simply supported and the unloaded 
edges (a) either both free or one edge clamped and the 
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remaining edge free. However, for plates with two loaded 
edges simply supported and the unloaded edges (b) both clamped, 
the buckling loads increase continuously from © = 0° to 0 = 

40° approximately whereafter the trend is reversed for four- 
and six “layered and orthotropic laminates. The variation for 
a two~layered laminate shows additional minima approximately 
at 9 = 20° and 70°. The difference in the response of 
buckling load between boundary conditions (a) and (b) can be 
attributed to the difference in the conditions on in-plane 
displacements. u° and v°, in addition to Poisson's effect 
which assumes predominance in the latter case. Another 
important observation for the type (b) boundary condition is 
that the buckling load at © = 90° is higher or lower than at 
© = 0°, This can be due to the occurrence of higher mode 
shape (i.e., m = 2 as compared to m = 1 at © = 0°) and/or 
the effect of ^'1/^2 noted that the buckling 

mode in the X“direction clianges from m = 1 to m = 2 at 
© i^;; 40°. This feature of the type (b) boundary conditions 
is distinctly different from type (a) boundary conditions 
wherein the buckling mode always retains the single wave 
shape, for the aspect ratios considered. 

The designer may also be interested in the optimal 
thickness ratio of the laminae leading to maximum buckling 
ibads. To this effect, the variations of buckling loads with 
thickness ratio for four- and six-layered laminates under 
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-various boundary conditions are studied. In all cases the 
buckling loads are found to be maximum at the optimal thick- 
ness ratios derived in Chapter 2 for various configurations 
of the laminates. However, it is of interest to note that 
the buckling loads are much more sensitive to thickness 
ratios for plates with clamped edges than for those with free 
edges. The use of optimal thickness ratios may, therefore, be 
more meaningful for plates with maximum restraints at the 
edges. The variations of buckling loads with cross~ply ratio, 
M for antisymmetric cross-ply plates with free edges, one 
edge free and one clamped and those with all edges simply 
supported are also considered. It is found that the increase 
in buckling load with increase in cross-ply ratio is maximum 
for plates with free sides than for plates with all edges 
simply supported. The maximum increase in buckling load is, 
however, limited to the cross-ply ratio range of M = 1 to 5 . 
Por a designer, this may indicate the optimal cross -ply ratios 
beyond which the increase in the ratio becomes unnecessary 
to achieve a given buckling load. 

An important aspect of the variations of buckling 
loads with plate aspect ratio for various boundary conditions 
is also considered. It is found that for plates with two 
loaded edges simply supported and the remaining edges free, 
the buckling load per unit width of the plate is fairly 
constant with aspect ratio and the buckling mode shape in the 
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x-dirsction always occurs at m = 1 for all aspect ratios. 

This is distinctly different from the other types of boundary 
conditions considered in this thesis. The cusps in these 
variations occur because of the change in mode shape and 
their location with respect to aspect ratio is found to be 
dependent upon (i) the number of layers (ii) the fibre orien~ 
tations and (iii) the boundary conditions. For a symmetric 
45° oriented plate, however, the location of cusps is found 
to be independent of number of layers. 

The buckling of composite plates under the action 
of both the in-plane and lateral loads and with two loaded 
edges simply supported and various other boundary conditions 
on the remaining edges is also investigated. The lateral 
loads are represented by the well-known Winkler-Pasternak 
medium. It is observed that both the material coupling and 
Winkler-Pasternak foundation stiffnesses effect the buckling 
response of the antisymmetric cross- and angle-ply plates 
significantly. However, the change in the Pasternak foundation 
stiffness exerts a far greater change in the buckling load 
than does a simi,lar change in the Winkler foundation stiffness. 
Furthermore, these effects are more significant for plates 
v/ith two loaded edges simply supported and the remaining 
edges free than for plates with other boundary conditions. 

The problem of the buckling bf antisymmetric cross- 
and angle-ply plates under the action of pure bending and . 
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with all edges siraply supported is solved employing the 
Galerkin techniq.ue. The variations of buckling load v/ith 
aspect ratio, number of layers, fibre orientation and thick" 
ness ratio of the laminae of a laminate ore studied » The 
general observations on the buckling response made with this 
type of loading are common with those of uniform uniaxial 
loads. However, the syrnmetry of buckling loads at 0 = 30° 
and 60° obtained for uniform uniaxial load v/ith all edges 
simply supported is absent for the corresponding case of pure 
bending, 

8 . 2 SCOPE HOR PURTHER WOBIC 

The optimal thickness ratios leading to special 
orthotropy of antisymmetric angle-ply laminates have been 
derived in this thesis. However, for general unsymmetric 
angle-ply laminates, partial reduction of coupling is shown 
to be possible with certain thickness ratios assuming an. 
arbitrary ply stacking s6q,uence. It may, therefore, be 
possible to work out the suitable combinations of ply stacking 
sequence and thiclmess ratios that may lead to complete 
elimination of coupling stiffnesses even for more general 
types of laminates. Such combinations will give additional 
choice of optimal thickness ratio and ply stacking sequence 
to the designer who may be interested in the optimization. 


201 


It is seen that for antisymmetric cross- and angle-" 
ply laminates considered in this thesis, it is possible to 
obtain the closed form solution to the various problems of 
buckling in view of coefficients Ag, A^ , A^ and Ag of the 
governing eighth order partial differential equation (3.1 4-) 
becoming zero. However, for more general types of laminates, 
this ajTialytical simplification is not possible and full 
stiffness matrices (A- B . . and B . . ) accompanying such ■ 
laminates can male© the analysis more complicated .• Pinlte 
difference or finite element techniques can be conveniently 
used to yield aai approximate numerical solution to such 
probl eras . 

Equation (3.22) can be used to study the buckling 
of laminates under combined normal and shear loads which is 
of direct practical interest to the designer. Interaction 
curves can be prepared from these results. Equation (3.22) can 
a].so be applied with suitable modifications to study the 
vibration of laminated plates with benefit. 

The study of the Influence of conditions on u° and 
v° on the buckling load variation with orientation in the 
case of free edges can also be of interest to reveal the 
relative importance of the various end conditions in deter- 
mining the buckling loads. 

' 'O: 
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